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pLn ' In the present article we are interested in further investigations for the barotropic 

compressible Navier-Stokes system endowed with a non-local capillarity we studied in 
[7] . Thanks to an accurate study of the associated linear system using a Lagrangian 
change of coordinates, we provide more precise energy estimates in terms of hybrid 
p ^ ' Besov spaces naturally depending on a threshold frequency ( which is determined in 

function of the physical parameter) distinguishing the low and the high regimes. It 
allows us in particular to prove the convergence of the solutions from the non-local 
, to the local Korteweg system. Another mathematical interest of this article is the 

' study of the effect of the Lagrangian change on the non-local capillary term. 



1 Introduction 

[ 1.1 Presentation of the system 

00 

• The local and non-local Korteweg systems aim to study the dynamics of a liquid-vapour 

, mixture in the diffuse interface approach (DI) , where the phase changes are seen through 

fSJ i the variations of the density. These systems are based upon the compressible Navier- 

' Stokes system with a Van der Waals state law for ideal fluids, and endowed with a 

capillary tensor modelling the behaviour at the interfaces between the phases. This 
capillary term was introduced in the DI approach in order to obtain physically relevant 
solutions by penalizing the high variations of the density. 
^ ■ We refer to [7] for a physical presentation of the diffuse interface model, and of the 

, local and non-local Korteweg systems. Let us recall that the local model of the capillary 

term was introduced by Korteweg and the non-local model was introduced by Van der 
Waals and renewed by F. Coquel, D. Diehl, C. Merkle and C. Rohde (for an in-depth 
presentation of the capillary models, we refer to [33] and [10]). 

Let p and u denote the density and the velocity of a compressible viscous fluid. As 
usual, /3 is a non-negative function and ti is a vector-valued function defined on W^. In 
the sequel we will denote by A the following diffusion operator 

Au = fiAu -|- (A -|- ^)Vdivu, with p > and v = X + 2fi > 0. 
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The Navier-Stokes equations for compressible fluids endowed with internal capillarity 
read: 

J dtp + div {pu) = 0, 

y dt{pu) + div {pu ®u) - Au + V{P{p)) = KpVD[p]. 

Let us mention that the capillary coefficient k may depend on p but in this article only 
the constant case is considered. In the local Korteweg system (NSK), the capillary term 
D[p] is given by (see [16]): 

D[p] = Ap, 

and, in the non-local Korteweg system (NSRW) (introduced in its modern form by C. 
Rohde in [32] and also [10], see Van der Waals [36] for the original works), if (j) is an 
interaction potential which satisfies the following conditions 

(|.| + |.|2)0(.) G i^(M'^), / (t){x)dx = l, (/> even, and 0, 

D[p] is a non-local term: 

D[p] =(j)*p- p. 

If we compute the Fourier transform of the capillary terms, we obtain {4>{C) — l)p(^) 
in the non-local model, and — [Cl^PlO in local model. 

We are interested in the closedness of the solutions of these models when (/>(^) is 
formally "close" to 1 — |^|^. In [7], we approximated the local Korteweg model {NSK) 
with a non-local model such as system {NSRW) where we chosed a specific function (j)^ 
in the capillarity tensor. In this paper we will once more consider the following non-local 
system: 

( dtPe + div {peUs) = 0, 
{NSRWe) j dt{peUe) + div {peU (g) Ue) - AUe + V{P{pe)) = Pe^V((/>e * Pe - Pe), 

where 

1 X 1 \^\^ 

'^^ = ^'^(7) ^ith cj){x) = -^e-— 

For a fixed ^ the Fourier transform of 0£ is (^£(^) = e"^^'^'^, and when e is small, 

is close to — 1^1^. 

We will consider a density which is close to an equilibrium state p and we will intro- 
duce the change of function p = p(l + q). For simplicity we take p = 1- The previous 
systems become: 



{K) 



dtq + u.Vq + (1 + g')div n = 0, 

dtu + u.Vu -Au + P'{l).Vq - nVAq = K{q).Vq - I{q)Au, 



and 

{RWe) 

( dtqe + Ue-Vqe + (1 + qe)divus = 0, 



K 



dtUe + Ue.VUe - AUe + P'{l).Vqe " {(j^e * Qe - Qe) = K{qe).Vqe - I{qe)AUe, 



£2 



2 



where K and / are real- valued functions defined on M given by: 

K(,)=(p'(l)-£^] and /(,)- ' 



l+q J q+l 



1.2 Existence results 



Let us now recall some results concerning the local and non-local Korteweg systems. 
As for the compressible Navier-Stokes system (we refer to [3], [11], [17], [6]) both of 
these systems have been studied in the context of the existence of global strong solutions 
with small initial data in critical spaces for the scaling of the equations. For example, 
concerning the strong solutions, we refer to [14], [22] (and [19] in the non isothermal case) 
for a study of (NSK) system, and to [18] for (NSRW). In [23], we show the existence 
of global strong solution with large initial data on the rotational part when we add a 
friction term. 

Let us mention that the well-posedness of the compressible Euler Korteweg system (when 
/I = A = 0) has been studied in the case of variable capillary coefficient by Benzoni, 
Danchin and Descombes in [4]. 

The solutions of the compressible Navier Stokes or non-local Korteweg systems have 
the same behaviour. Namely the density regularity is separated by a frequency threshold: 
in low frequencies, the solution is subject to a heat-type smoothing, and in the high 
frequencies, there is only a damping effect due to the term of pressure (modulo that 
the pressure is at least locally increasing with respect to the density). The solution 
of the local Korteweg system is more regular: for all frequencies, we have a parabolic 
regularization on the density (see [14, 22]). We refer to the appendix for the definitions 
of the Besov spaces introduced in the following results. 

def 

Theorem 1 ([14]) Assume that P'{1) > 0, min(;U,z/) > where v = 2/x + A, that the 
initial density fluctuation qq belongs to B21 H B21, and that the initial velocity uq is in 

{B21 )'^. Then there exist constants rjK > and C > depending on k, h, v, P'(1) and 
d such that if: 

then system (K) has a unique global solution {p, u) such that the density fluctuation and 
the velocity satisfy: 

q G C(M+, bI~^ n b|i) n Li(M+, bI^^ n bI^^), 

u G C(M+,s|~VnLi(M+,^|^V- 

Moreover the norm of {q,u) in this space is estimated by the initial norm C{\\qo\\ d_-^ d + 

Ikoll ._,). 
^2,1 

Remark 1 Further in this article R. Danchin and B. Desjardins provide a Fourier study 
of the linearized system and observe different behaviours whether the quantity — 4k 
is positive, negative of zero. In all cases they obtain parabolic regularization. 
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Remark 2 In [21], the second author obtains some generahzations of [14] in as much as 
with a specific choice on the capillarity and the viscosity k{p) = ^, fi{p) = p, we obtain 

the existence of global strong solution with uq E -62^00 ^^(^ + ^o) £ -^2^00 ^ ^2oa 
(let us point out that in this case we can work with discontinuous initial density). We 
also refer to [23] for the existence of global strong solution with large initial data on the 
irrotational part when we added a friction term. In the sequel as we shall deal with a 
Lagrangian change of coordinate which requires a Lipschitz control on the velocity, we 
need to work in the framework of the functional space introduced in [14]. The second 
reason is that the existence of global strong solution for the system (NSRW) with initial 

d d -j^ 

data in Bi^^ n -Bl 00 ^^^^cLins open and is probably false in general. 

As for the compressible Navier-Stokes system, in the (NSRW) model the density fluc- 
tuation has two distinct behaviours in some low and high frequencies, separated by a 
frequency threshold. This naturally leads to the definition of the hybrid Besov spaces, 
involving two different regularities for low and high frequencies, introduced in [7] and 
defined for = ^og2{-^j^) — 1] (7 is a constant and Co = |) and s,t E M by : 



j;2'^||A,g||i.+J^l2^iAzg|U2. (1.1) 



I II '^^Z 
g Ds,t — 

1<U l>k 



Definition 1 ([7]) The space is the set of functions (q, u) in 
(a(M+, n Bl,) n L\R+,BI+''' n Bl+^'')) x (Ch{R+,B^') n l\r+, 501))'' 



endowed with the norm 



de f 

\\{Q,u)\\ei = lhllL<x>S--i + IklliocB;'-! + lkllL°°B, 



2,1 



We first state the global well-posedness for system (RWs) with uniform estimates with 
respect to e (see [7]): 

Theorem 2 ([7]) Let e > and assume that min(/i, 1/ = + A) > 0. There exist 
two positive constants i]r and C only depending on d, n, p., A and -P'(l) such that if 

go £ B2 ^ n B2,\j "Uq G -82^,1 and 

Ikoll . #-1 . # + Ikoll . d-i < TlR 

d 

then system [RWe) has a unique global solution (psjUe) with {qe,Ue) E Ei such that: 

\\{qe,Ue)\\ d <C{\\qo\\ d +||no||.d J. 

nBl, Bl, 
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Remark 3 Note that in the low frequency regime (I < le ^ — ^logej, the parabolic 
regularization for is the same as for the Korteweg system, indeed the low frequencies 

ofqe are in n ^|^^). 

The main result in [7] is the fohowing: when the initial data are small enough (so that 
we have global solutions for (K) and (RWe)) the solution of (RWe) goes to the solution 
of (K) when e goes to zero. 

Theorem 3 ([7]) Assume that min{ij,,2ij, + X) > 0, P'(l) > and that qo e nBl^, 
no G B21 ■ There exists < r] < min(?7x, 77/^) such that if 

Ikoll .4-1 d + \\uo\\ . d_i < 77, 

then systems (K) and (RWe) both have global solutions and \\{qs — q,Us — u)\\ d tends 

to zero as e goes to zero. Moreover, with the same notations as before, there exists a 
constant C = C{ri, k, P' {!)) > such that for all a g]0,1[ (if d = 2) or a g]0, 1] (if 
d > 3), and for aU t G M+, 

\\{q,-q,Ue-u)\\ <C£", 



This results relies on the following estimates: 

Proposition 1 ([7], Proposition 1) Let e > 0, s G M, / = [0, T[ or [0,+oo[ and v G 

L^{I,B^^ ) n L2 (/, ^1 1 ) . Assume that {q,u) is a solution of System (LR^) (see below) 
defined on L There exists a constant C > depending on d, s, fi, u, p, n such that for 
ah t G /, 

TooR^^l ^ 7oo rs-1 + K7 Too RS + 71 ~I~ 9 TIR^ + I''^ ^ 7l R'' + 2'^ 

C/o{liV»;(r)|| , +|b(r)||2 , )dT 

+ K 71 R»-l + 71 Rs + 7lRS-l)- 

^t^2,l ^t^2,l ^t^2,l ) 



Remark 4 Let us point out that the main difficulty consists in obtaining the previous 
accurate estimates in Besov space depending on the parameter e for the linear system 
(LRf,). Indeed it is quite tricky to deal with the convection terms (let us point out that 
it is absolutely necessary to integrate the convection term v ■ Vq in (LRs) in order not to 
loose regularity on the density in the remainder term F, indeed (RWs) does not provide 
any regularizing effects on the density in high frequencies), in [7] we use energy methods 
and symmetrizers. 

The goal of this paper is to propose a more robust method which allows us to precisely 
keep track of the dependance with respect to the physical coefficients: viscosity and 
capillarity. 
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1.3 Statement of the results 



Classically, as in the study of compressible Navier-Stokes systems-type in critical spaces 
(see [11, 6, 17]), for proving the theorems 2 and 3 (see [7] section 2) the key point consists 
in obtaining a priori estimates on the following advected linear system (e > is fixed 
and for more simplicity we write {q,u) instead of {qs,Us)). 



{LRe) 

With 



dtq + v.Vq + div u = F, 



dtu + v.Vu — Au + pVq oV{(t>e * q — q) = G. 



Au = fiAu + (A + ^)Vdivit. 

As mentioned before, in the present article we are interested in obtaining, from a 
different point of view, these energy estimates by using a Lagrangian change of coordinate. 
Our method provides a more precise dependency of the various constants with respect 
to p = P'{1) and the well-known ratio v'^ /Ak (that also appears in the local system, and 
in any evanescent viscosity-capillarity limit. In the ID case, we refer to [8] for the study 

2 

of this limit (when = 0(1)). This result is proven via the introduction of an effective 
velocity) . 

Remark 5 In Proposition 7 we provide two equivalent (and more accurate) expressions 
for the hybrid norm ||.||^s,s+2 and we will only use them: 

II/IIbi-'^ - II " b| , - E ^ " 'J ' ^ 2^11 A,/|U2 

~5]min(l,22j)2^-^||A,/||i2 (1.4) 

Consequently the non-local capillary term j-^g^g [j^ fg.Qi same regularity as 

the local capillary term VAq : both of them belong to L\{B2i H i?! i )• 

Let us now give the main result of the present article, which is a sharper version of 
Proposition 1 (see the previous section): 

Theorem 4 Let e > -^ + 1 < s < ^ + 1, I = [{),T[ or [0,-Foo[ and v G L^{I,B^^^)r\ 
■ i 

L^{I,B2i)- Assume that {q,u) is a solution of System {LR^) defined on I. There exists 
eo > 0, a constant C > depending on d, s such that if e < Eq, for all t & I (denoting 
u = fi + 2X and uq = min(i/, /i) ), 

II"IIl-b»-i + Ikllz-B^-i + ^h\\LfBi^ +MML]Bif + ^lkllziB|+i-''-i + i'^lkllzii?i+2>'' 

C 2C„,, A||V^;(r)|| . +\\v{t)\\\, )dT 

X (ll^^ob-1 + lkolU-1 +j^lkobs_^ + II^IIljs-i +i^ll^llzji?i_, + ll^llzjB-0- ^^-^^ 
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where 

1 4k u"^ 

= '^max(VP,— )max(^,( — 

_ l + |A + ^|+/i + z. 1 

^msc — r nidXI i, — ^ I 

Remark 6 The viscous coefEcient Cmsc satisfies: 

C l+2i' 



i±^+max(l,^) JfA + /i>0, 
i±^+max(l,i,) JfA + ^<0 



and when both viscosities are small, we simply have Cmsc < max(l, ■\). 



Remark 7 The norm \\q\\ri ns+i.s-i in the left-hand side has to be compared to the term 
llgll^-igs+i.s from proposition 1. Thanks to the last term in the left-hand side, there is no 



loss of regularity. We only chose to write it this way in reference to the more meaningful 
equivalent expression from (3.49). 

Remark 8 Let us give a few comments about the advantages of the Lagrangian method 
compared to the symetrizers techniques used in [7]: 

• First we are able with this method to provide accurate estimates tracking the 

2 

physical coefficients and especially the influence of the ratio |^ (we also refer to 
[14] for the importance of this ratio in the local Korteweg system). Indeed it plays 
an important role when considering the vanishing viscosity-capillarity process. In 
the one dimensional case we proved in [8] the global convergence of the classical 
Korteweg solutions to the global weak entropy solutions of the compressible Euler 
system when k = and v goes to zero. Let us also mention the works of Lax and 
Levermore ([30]) who consider the non-viscous case and prove a vanishing capillarity 
process from the KdV equation to the Burger equation in the context of dispersive 
shock solutions (the main tool is the inverse scattering theory). Letioch also studies 
the KdV equation with viscosity, and shows that the previous ratio is critical for 
the convergence towards a weak entropy solution or to a dispersive shock solution 
(notice that this also seems to be observed in numerical simulations). In particular 
we expect our linear a priori estimates to be useful to study the vanishing process. 
This is the object of a future work. 

• Seconds, the coefEcients of the last two terms in the left-hand side of (1.5) are v 
and , which in the setting k ~ i/^ with v small, are larger than the obtained 
in [7] (see the end of section 2.1 therein). 

• Another important feature of the present paper is that this method allows us to 
deal with Besov spaces constructed on general L^ spaces with r ^ 2 in the spirit 
of [6]. Indeed, such results cannot be obtained with symmetrizers methods, which 
are by nature based on energy methods and scalar products in L^. We refer to the 
appendix for results in this direction and to [7] for the compressible Navier-Stokes 
system. 
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The paper is structured in the fohowing way: instead of using energy methods and 
symmetrizers (see [11, 7]), we wih first consider the hnear system (LR^) without transport 
(this is the object of the second section, where we precisely study the different frequency 
regimes). In section 3 after presenting the Langrangian change of coordinates, we provide 
equivalent expressions for the hybrid Besov norm, and perform the Lagrangian change 
of variable, as introduced by T. Hmidi in [24] (in high frequency regime, transport terms 
prevent any direct use of the linear estimates). 

As in [15, 6, 25, 26, 27] we get estimates on the advected linear system and for this we 
need to bound additional external force terms generated by the change of variable. This 
is done thanks to estimates dealing with the action of a lagrangian change of variables 
on frequency truncation, such as the one proved by Vishik (see [37]), we also refer to 
[3]. In this part we focus on the main difficulty of the present article which consists 
in estimating the commutator of our non-local capillary operator under the Lagrangian 
change of variables. We end this section by giving an extension of our estimates allowing 
results in Besov spaces defined on L'^ with r ^ 2. The appendix is devoted to give the 
main tools of the Littlewood-Paley theory and recall classical results on the Lagrangian 
flow. 

Remark 9 Let us mention that another approach would consist in introducing an effec- 
tive velocity as in [1 7] in order to diagonalize the system in a certain way and to cancel 
out the coupling between density and velocity. This is the object of another article. 

Remark 10 For general considerations about non-local operators we can refer to the 
work of Rohde and Yong (see [34]) and the recent paper of Alibaud, Cifani and Jakobsen 
(see [1]). 



2 Linear estimates 



The aim of this section is to obtain linear estimates for the following system: 

{dtq + divn = F, 
dtu - Au + pVq 9 V(0£ * q - q) = G. 



(Le) 

With 



£2 



Au = jsAu + (A + ^)Vdivu. 

Let us state the frequency-localized result that we will use in this article: 

Proposition 2 Let e > 0, s G M, / = [0, r[ or [0, -|-oo[. Assume that (g, u) is a solution 
of System (L^) defined on I. There exists Eq > 0, a constant C > depending on d, s, 
Co and Co such that if e < Eq, for all t £ I (as usual uq = min(i/, ^u) ), and for all j G Z, 

||Aju||ioo^2 + fo22j||Aju||ii^2 + (1 + 1^2^) (^\\Ajq\\L^L2 + umin{^,2'^^)\\Ajq\\^L2^ 

< Cmax(vp, ^)max(^, (^)^) 

X ((1 + v2^)\\AM\l^ + l|Ajno||L2 + (1 + v2:^)\\kjF\\^r^2 + \\AjG\\i^.^^ (2.6) 
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2.1 Study of the eigenvalues 

As in [11] or [6] we first introduce the Helmholtz decomposition of u. If the pseudo- 
differential operator A is defined by A/ = /(.)), we set: 

{V = A~^divti, 
_i (2-7) 
w = A curlu 

then u = — A~^Vf + A~^divu' and the system turns into: 
dtq + Av = F, 

(L^) \ dtv - uAv - pAq + ^A{(t)e * q - q) = A"Miv G, 

dtw — fiAw = A~"'^curlG. 

The last equation is a heat equation, easily dealt thanks to classical heat estimates in 
Besov spaces (we refer to [3] chapter 2), so we can focus on the first two lines and compute 
the eigenvalues and eigenvectors of the matrix associated to the Fourier transform of this 
new system: 

dtq+m=F, 

As the external forces appear through homogeneous pseudo-differential operators of de- 
gree zero, we can compute the estimates in the case F = G = and deduce the general 
case from the Duhamel formula. So we finally study the following system: 

dt( l)=Ai0( I) with AiO:-^ ' 



The discriminant of the characteristic polynomial of the matrix is: 
A = |eP (.^|eP-4(p + ^(l-e-^l«l^) 

and thanks to the variations of function : x i— )■ v'^x — 4 + ^(1 — e~^^^)^ , we obtain 
the existence of a unique threshold > such that 

f<oif \e<xe, 

^^^\>oif iep>x,. 

2 2 

We emphasize that when > 1, is an increasing function on and when |^ < 1, 
is decreasing in [0, log(^)] and then increasing. 

Proposition 3 Under the same assumptions, we have: 



where a^^u = a(^) is the unique positive root of function x i— t- 1 — Moreover 
we have ^ - 1 < a«;,i. < ^. 

Proof: First, as the function h : x ^ ^— | — decreases from [0, oo[ to ]0, 1], we easily 
prove that the following function: 

=^— = 1-^ 2 2 

is an increasing function from ]0, oo[ to ] — oo, 1[, and for a fixed a; > 0, it is increasing 
with respect to e (Vx > 0, VO < e < e', 5e(x) < gi£'(x)), which implies that e i— Xg 
increases when e decreases to zero. 

Then, computing ge{-[;^) = e — which is positive when e is small enough, as 
is increasing, we obtain that x^ < when e is small enough. Computing (?e(^) = 
— ^/i(e2^) < we easily get that Xg > Finally we obtain that: 



• If 1^ > 1, as for all x, < h{x) < 1, we have 1 - ^ - ^ < Qeixe) = < 1 - ^ 
which implies that the sequence x^ is bounded (and therefore convergent because 
monotonous), so that e'^Xe goes to zero, and thanks to the fact that /^(xe) = and 
h(x) — )• 1, we obtain that x^ — )■ ff. . 

• If 1^ = 1, guided by the value of ^^(y) = 1 — — h{£C), and thanks to function 
study and Taylor expansions, we get that: 

A/ < Xe < ' 

V K e 

• If 1^ < 1, guided by the value of fl'e(p-), we obtain that for sufficently small e: 

where a(f^) is the unique positive root of function x i— )■ 1 — f — and satisfies 

Then we can compute the expressions of g and v: 

-For the low frequencies (A < 0), when |,^| < ^/x^, we have: 

m = I ((1 + sw)e''- + (1 - ^)e*^-) fo(e) - ^^tei^-'oco, 

with: 

5(0 = ^TMm = + (2.9) 



and 



z^2|^|2 V" e 
A± = -^(l±i5(e)). 
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-For the high frequencies (A > 0), when |,^| > -y/xj, we have: 

m = ^ ((1 - T^y)^*'" + (1 + i^y)^*'-) *'(^) + '^^M^^'^'(^)' 



with: 



^(C) = = - ^ (P + ^(1 - (2.10) 



and 

A± = -^(l±i?(0)- 



Remark 11 ft wiJ] be crucial for the time integration to observe that 



p + 5(1 - ^"^''^'') = ^(1 - ^(0)(i + m))- 



Remark 12 Jii the low frequency regime, both eigenvalues provide parabolic heat regu- 
larization. In the high frequency regime, when |^| is large, A+ ~ — z^lCP (which generates 
parabolic regularization) , but A_ ~ — (p+ t-), that only provides a damping. This is the 
same behaviour as for the compressible Navier-Stokes system, and we refer to [11], [3] 
and [6] (for a precise computation of the Fourier transform of the solutions). 

Remark 13 These results have to be compared to the same study for the local Korteweg 
system, in this case the matrix becomes: 

-lei 

And the discriminant of the characteristic polynomial, A = ((i/^ — 4k)|^P — Ap^, 
satisfies: 

. If g < 1, A < and A± = -H^P^^I^I V4p-(^^^ 4k^ (parabolic regularization 
everywhere), 

• 7f|^>l,A<0 then A > with threshold ^^^^^ , and in the high frequency regime, 

when 1^1 is large, A_ ~ — ^^^-(1 — — ^) which, contrarily to the previous cases, 
generates parabolic regularization. 
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2.2 Pointwise estimates 
2.2.1 Thresholds 

As seen previously, we recall that satisfies: 

( 4p 



i,/^ if^ = l 



«(g) 

£2 



But in the following we will be able to get a frequency threshold of size l/e^ in each of 
the three cases, for that we will not directly use as the announced threshold but a 
larger term, y^, defined the following way (we refer to (2.8) for the definition of Qs): 



where M = —. 

ye \ 2M' - 4' 



(2.11) 



As Qs is an increasing function from ]0,oo[ to ] — oo, 1[, we have < y^. The following 
property shows that i/s is large enough: 

Proposition 4 Under the previous estimates, there exist two constants < 71 < 72 and 
eo such that for aU < e < Eq, 

Xs < ^ <ye < ^■ 

Moreover, if M = ^ > j, then 71 and 71 are universal constants, and if M = |^ < | 
then we have: 

11 3 

- < 1 < a(M) < 71 < 72 < — . 



Proof: Here it will be more efficient to compare f;^ to | instead of 1: 
First case: If M = |^ > | we define 71 and 72 by: 

^(71) = ^ and /i(72) = ^. 
As gei^) = 1 - - jih{C), we obtain that 



'21) = 1 _ J 4^,2 ^i__L 

'e'^' 2M 1/271 2M 



and if e is small enough (depending on p, u, 72 and M), 
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Second case: If M = ^ < |, 71 and 72 are defined by: 

3 1 
/i(7i) = -M and h{j2) = ^M. 

Like previously, we compute: 

5'^ 72 = 1 ~ 7 ~ ^"'^ < 7 = fi'e U/e), 
4 71 4 

and if e is small enough (depending on p, ;/ and 72), 

S-e — = 1 - - ^— e > 9e{ye), 

2 z^^72 

In both cases we conclude using that function is increasing. 

In the second case, to get the lower estimate we use that h{-yi) < M = h{a{M)) to 
obtain that (remind that the function h is decreasing) a(M) < 71. We conclude thanks 
to the fact that a(M) > ^ - 1 > |. ■ 

2.2.2 Estimates 

Now that we have defined the frequency thresholds and y,. we can state the main 
result of this section: 

Proposition 5 Under the previous notations, there exists a constant C, such that for 
all j and all ^ G 2^C where C is the annulus G M'',co = | < |^| < Cq = §}, we 
have the following estimates (we denote by fj = Ajf and we refer to the appendix for 
details on the Littlewood-Paley theory): 

. If\(\<^: 

2 '2j 

(1 + u^MjiOl < Ce-^ ((1 + u2^)\q^,{0\ + (1 + ^)l^^(e)l) , 

. m)\ < Ce-^ ((1 + ul^il + VP)(1 + ^)|5^(0I + 1^(01) • 
• If ^ < 1^1 < 

' (1 + y^mm < T^e-^(i-™) ((1 + ^2^)1^(01 + (1 + > 

^v,{0\ < T^e-^(i— ) ((^2^-)|5^(OI + 1^^(01) , 
wLere m = i if M = g < |, m = if M > |. 

. If\C\>^: 

(1 + u2^)mO\ < C (e-"^ + e-^(^-^"'^)*) ((1 + z^2^)|g^(e)| + (1 + -^)\^^m) : 

< C (e-^ + (1 - 7^;(^)e-^(i"V^^(^o2^^))) (^2i|^(e)| + |5^(C)|) . 
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Proof: We split the study according to the two frequency thresholds: 

Low frequencies: assume that ^ G 2^C with |^| < ^/x^, then the discriminant A < 
and the system has two non-real conjugated eigenvalues, thanks to (2.9), we can write: 

with: 



and 



A± = -^(l±i5(0). 
Here we have to cope with two difficulties: 

1. neutralize S{S,)~^ for frequencies close to -y/xj, as 5'(^) goes to zero as |^| goes to 



2. if we are not careful, the first term in the expression of the velocity will provide 
either e~^, either 2"^-' which would make our estimates useless. 

The first point will be adressed by considering the following block: 



-e 2 iut\£\' = -iiyt\£\'e 2 



As for all x > 0, xe~^ < -e~^/^, we obtain that: 



gt\+ _ gtA- 4 ^tie|2 

We then deduce: 

\m)\ < Ce-"^ {\q^M)\ + ^l^(^)l) • (2-13) 

And finally: 

y2%{£,)\ < Ce-^ + 1^^(01) ■ (2.14) 

For the velocity, we can also use (2.12), except for the first term (that provides large 
coefficients described as difficulty 2 in the previous page.): 

A = I [p+—{l - e"^ 1^1 



We can observe that S{£) > y i;^|p ~ ^ ~ u]^ Y^4p — i^^l^p, which makes sense only 
when 1^1 < l^Jp/v. Therefore we split the study into two cases: 
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If 1^1 < then v\i\S{i) > v^4p-z^2|^|2 > ^ ^nd as 1 - e"^ < x, we get: 

V2p V2p 



Using this together with (2.12), we get: 

going back to the density we can estimate the last term, and in addition to (2.14), 
we obtain: 



mo\<ce — ^(15^(01 + 1^^(01)- 

If ^/2p/u < ICI < y/x^, then we have gei^) < gs{\C\'^) < Oeixe), which imphes: 



4 / K , =.2|f|2,\ , 4k, 



and we obtain for the velocity: 

|^;.(0| < Ce-^ ((1 + ^)u2^q^{0\ + \v^m) ■ 
and for the density fluctuation, thanks to (2.13), as |^| > yflpjv^ we also have: 

\qAi)\ < ce — — + ^\v^m) ■ 



Finally, gathering these estimates gives the first point of the proposition 5. 
Remark 14 Note that ge(^) < dei^) < = Qeixe) so ^ < in any case. 

Remark 15 We emphasize that in the case of the compressible Navier-Stokes system 
(see [6]), such a distinction was not necessary as the analogous of 5(0 was y'^^^p^^^ 
and the threshold was located at = 

High frequencies: assume that ^ € 2^C with |^| > > ^/x^, then the discrim- 
inant A > and the system has two real eigenvalues. Thanks to (2.10) we write the 
localized quantities: 

\m = - (p + ^(1 - "^iy^^&O + ^ ((1 + T^y)^*'" + (1 - Tm^^'"-) ^ 

with: 
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and 

A± = -^(l±i2(0)- 

We have e + < e 2 and as explained, when |^| is large we cannot hope to get 
parabolic regularization from A_ . But the definition of implies that in this case 



9em^) > 9e{ye) = { 



, if M = 

l-2M>|ifM>|, 



therefore i < R{^) < 1, 1 + R{0 G [1,2] and 1 - S [0, i]. More precisely, as 

I^P > ?/e > 77, we have: 



So that we can bound A_ = —^^^{1 — R{£,)) and e^^^ < e ^ which allows to 

estimate the density fluctuation the following way: 

< C{e-'-^ + e-'*<i-=""') (|5SJ(£)| + ;;^|SC(Ol) ■ (2.16) 



Using successively that |^| > co2^ and |^| > ^/x^ > yields the estimate given in 
proposition 5. 

Concerning the velocity, we have two main difficulties in this frequency domain: 

• Obtain parabolic regularization for the velocity, even when the "damping eigen- 
value" A_ is involved. 

• As previously, estimating like before the first term provides coefficients (un- 
bounded in e) or 2^-' (too many derivatives). 



In order to neutralize these large coefficients, we rewrite the velocity as in Remark 11: 

j-A = (1 - i?(e))(i + Rior^^j^iq^ioi 

1^.(^)1 < A + IB + C with J]B= i(l + -^)e*^+|5^(0|, 



2(74) -1)^*'- 1^(01- 



We can easily estimate B: 

B<Ce "—\v^{0\- (2.17) 

The other two terms involve A_ and have to be handled carefully as we wish to get 
in time estimates: 

c < c(i - ii(e))e-^(^-^(«»i5^(e)|. 
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Using the fact that function is increasing, as cqP < |^| < Co2^ , we have: 



and 



1 - JgsiCp^^) < 1 - R{0 < 1 - Jgeicp'^), 



so that we can write: 



c < c (i - y^^^) 1553(01 



(2.18) 



We estimate the last term by: 



A < Ai + A2 with < 



A2 is easily dealt: 



" ^ ' 



Ai = ^(l-i?(0)(l + i?(6)- 



" e 2 



A2 = ^(l-i2(0)(l + i2(C))- 



^(01, 



A2 <Cz.2^e ^15^(01, 



and finally we get: 

A, < ( 1 - y^^^) e-^{-^5^)|S3K)|. 



Gathering these estimates together with estimates (2.17) and (2.18) leads to the an- 
nounced result in this case. 



Remark 16 Though small, terms like 1 — \J gsiCp"^^) will allow to get in time 
estimates with parabolic regularization. 

Medium frequencies: assume that ^ £ 2^C with < < We use the 

same formula as in the high frequency case. Due to the fact that R{(,) goes to zero when 
1^1 goes to we write the localized density and velocity as in the low frequency case: 



v,{0 = !^{l-Rmi + RiO) 



2 \ ^ - ' R{0 



■^^(6> 



RiO ' 2 

In this frequency domain, we define m such that: 

= y"^) < RiO = V^MF) < m =^ 



mm) 



1 / pt^- 



^(0 



i if M = — < - 



l__i^ifM>|, 
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As m < 1 we have 

[l-m<l- R{^) <l 

The aim is to get parabohc regularization and neutrahze the possibly vanishing -R(^) 
near ^/x^, by considering the following blocks: 



e 



As A± = -^(1 ± (see 2.10), we have: A+ - A_ = -i^\S,\^R{(,) and then (thanks 

to the variations of function h and to (2.19)): 

Finally, using once again that xe~^ < |e~^/^, we obtain: 

< — :— e ^(1-™). (2.20) 



' R{0 l-m 
From this we easily get that: 

C -l^fin ^ / 1 \ 

1^(01 < Y^^'^^ ^""^ + I^I^(^)lJ • (2.21) 

Using successively that ^ G 2^C and |^| > > the estimate of proposition 5. As 
A+ > A_, with the same method, we obtain the corresponding estimate for the velocity. 
This ends the proof of the proposition. ■ 

2.3 Time estimates 

As in the case of the compressible Navier-Stokes system (see [6] section 3.1), due to the 
choice Co = 3/4 and Cq = 8/3 (see the appendix), we can observe that there exist at 
most two indices = Jo " 1 or = Jq such that ^ e 2J[co, Cq] for j G {j^, Jq}. 
The aim of this part is to prove the following result, which implies Proposition 2. 

Proposition 6 Under the same assumptions as in Proposition 2, there exists a constant 
C such that for all j £Z (denoting M = ^): 

• For aU j < Jq, 

lbillL?°L2 + l^22^lbillLiL2 + (1 + 1^2^) (lkjllL-L2 + l^22j|kjllLiL2) < 
Cmax(^,M2) (^(1 + i.2^-)(l + ^)\\qo,jh^ + (1 + -^)\\vo,j\\l^ 

(2.22) 
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For all j > Jq, 



\v 



1 

(2.23) 



C7max(l,M) ( (1 + z.2^-)||goj 11^2 + (1 + — )||7;oj Hl^ 



Remark 17 Notice that due to the separated results in low and high frequencies, the 
coefficients here are more precise than in Theorem 4. 

Proof: As previously, three cases have to be considered: 

High frequencies: if j > j'q + 1 then ^ 2-' [cq, Cq] and for all ^ G 2^C, we have 
^ < y/y^ < Co 2-' < 1^1 < Co 2-' so that we only use the high frequency case from the 
previous proposition. Integrating with respect to ^ (we recall that the frequencies are 
localized in 2^C), we obtain thanks to the Plancherel formula: 

'(1 + i^2^)\\qj\\L2 < C (e-^ + e-^(i-^"'^)*) ((1 + iy2^)\\qojL^ + (1 + -^)\\voM^) ^ 

^ Ml^ < C + (1 - V^;(^)e-^(^-V^^(^ (-2i,o.||.. + Ml^) 

This immediately implies that: 

(1 + vT)\\q,\\LrL^ < C ((1 + y2i)\\qoAL^ + (1 + . 
hjWb^L^ < C (l/2^||goj||L2 + ||?;o,illL2) . 

The L|-estimates require a little more work: 

(1 + vV)h,\\^,^, S C (|£ + ^ri^TTi) ((1 + i'2i)\\m.j\\L' + (1 + 75)ll"0jlll») . 

""ill..'.' s (t;^^ + :*i^f^P) + . 

Thanks to the definition of 71 (see Proposition 4) we can write: 



so that: 



71 |2MifM = |:>f, 

2 



l-e~7. l|ifM = £<|, 



2e2 i.e2 ^2 / ^2 X ^2 ^ 

— + ^:^ 3^ = — 2 + - % <_(2+ ). 

1^71 k;(1 — e 'ij 1/71 \ K 1 — e '1 / 1^71 4 

Proposition 4 also implies that I/71 is bounded from above by a universal constant, so 
that we obtain: 

(1 + v2i)\\qj\\LlL^ < max(l,M) (^{1 + iy2^)\\qo,j\\L2 + (1 + -^)\\voJ\l^^ ■ (2.24) 
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We now turn to the velocity: as function in increasing, we get that 5e(cg22j) < 
5£(Cq2'^-') and then ) ^ ^ This term has to be bounded uniformly in e and 

j. Thanks to the expression and variations of we first write: 



X 



p 1 — e 



-£ c; 



22j 



r<2 — I 5 r'2 

<4ro_K ei <4ro /2 25^ 



as function x i— ?• 1 — e ^ is increasing. This yields: 

WvjWb^L^ +^'^'^^\\vj\\L\L^ < C {u2^\qoj\\L2 + ||t;o,j||i2) , 
so that estimate (2.23) immediately follows. 

low frequencies: when j < j^ — l we know that ^ [co2^ , Co2^] and here we will 
have to consider both cases |.^| < and < |^| < Ue- 



hjh^ < c (^||^;-1{|5|<^}||l2 + \\QA{^<\i\<^}\\L^ ) ■ 

Using the medium and low frequencies estimates from Proposition 5 provides: 
(1 + 1^2^)1^,11^2 < C i^e-^^ + ^— -e-^(i— )J 

X (^(1 + u2^)\\qojL^ + (1 + -^)|boj||L2) , (2.26) 

which implies: 

(1 + iy2^) {\\qj\\L^L^ + i^2^'\\qj\\LlL^) < 

^^^^ ((1 + i^2^)\\qo,j\\L^ + (1 + ^)lbo,,||L2) . (2.27) 

Let us recall that from Proposition 4 we have: 

1 (2 \f M — !^ < ^ 

1 < —— < I 4k - 4' < 4niax(l,M). 

~ 1-m |^4M if M > I 

The same can be done to the velocity and gives: 

Cmax(l, M2) (^(1 + z.2^)(l + + ^)\q^m + I^COl) • (2-28) 
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Which imphes 2.22 in the case i < Jq — 1- 

Threshold frequencies: in the case j G {jQiioli know that y/y^ G 2-'[co,Co] 
which may also be true for Therefore we are forced to write (even if the first term 

of the right-hand side is zero when M > 1...) 



hjh^ < C (^||g5-l||5|<^}||i2 + ||g5-l|^<|^|<^}||i2 + ||g5-l||5|>^}||i2 
Using all three estimates from Proposition 5, together with the specific bounds 



2L < £222, 72 



7l_ . .2n2,- , 72 
-12 

we end up with: 



(l + i/2^)||gJ|i2 <C e ^(i-'") +e"^^^'o^ M 

\ 1 — m / 

X (^(1 + i^2^')|koj||L2 + (1 + ^)lko,,llL2) , (2.29) 

The rest of the proof follows the lines of the low frequencies case, except that for the 
time integral, we get: 

i'2='(l+i'2')||«j||i,is < 

"((T^ + I^r^) ((i + "2')ll».ll^' + (i + ;;|)l"«.l^=). 

From Proposition 4, we have: 

72 < [2^ if M > f, ^^^^ ^ ^ fc if M > f, 
1-^-"^ ~\mf,^iM<l 7i-\T^<12ifM<|. 



which allows to write: 

1 _^ 72 



(1 - ?n)2 4ft; 1 - e" 



■71 



< C ^max(l,M2) + Mmax(l,-^)^ < Cmax(l,M2). 



Doing the same frequency truncations for the velocity achieves the proof of the proposi- 
tion. ■. 

2.4 End of the proof of Proposition 2 

Let us recall that in 2.7, w satisfies a classical heat equation, then for all j G Z (still in 
the homogenous assumption): 

Wwjh^L^ + ^^'i'^^\wj\\LlL2 < C||u70j||l2, 

and, together with Proposition 6 (we recall that = min(z^, n) = min(/i + 2A, fi) > 0), 
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• For all j < Jo, 

Cmax(^,M2) (^(1 + u2^){l + ^)\\qojL^ + (1 + -^)\\uoJl^^ , 

(2.31) 

• For all j > Jo, 

\\uj\\L^L^+uo2^^\\uj\\LiL^ + (1 + 1^2^) (^Ml^l^ + umm{^,2^^)\\qj\\LiL2^ < 
Cmax(l,M) (^(1 + u2^)\\qojL^ + (1 + -^)||noj ||l2^ , 

(2.32) 

which implies Proposition 2. ■ 

3 Proof of the advected linear estimates 
3.1 Presentation of the difficulties 

The aim of this section is to prove Theorem 4. A natural idea is to use Proposition 2 
and put the advection terms as external forces: we can write that for all j G Z, 

\\Aju\\L^L^+uo2^^\Aju\\^L^ + {l + u2^) (\\Ajq\\L^L^+umm{^,2^^)\\A,q\\LiL2) < 



C ^2 



(1 + ,.2^)\\A,qo\\L2 + ||A,i;o||L2 + (1 + ^2^) ( l|A,F||^j^2 + \\A,{v.Vq)\\^.^, 



+ (\\A,GhiL^. + \\A,{v.Vu)\\LiL^ 



. (3.33) 



The next step consists in multiplying by 2^^^~^') with — ^ + I<s<^ + 1 and sum over 
all frequencies. There are three additional terms to estimate. Thanks to paraproduct 
and remainder laws (we refer to (4.125) in the appendix) we have: 

||v.Vu||^.--i < ||T;Vn||^.-i + ||Tv«v|lij»-i + \\Riv,Vu)\\^s-i 
<c[\\v\\l^\\Vu\\^.-^ + \\Vu\\ d\\v\\ d +11^11^0,^^ II Vn||^.-i <^ll^ll^| 

\ J^oo.oo ^2,1 / ^2,1 

ill K 

<C||?;||.d ||ti||L-i||^i||L+i < TTT^i^oll^^lU^+i + C^t;— ||^'|P d IklU'^-i, (3.34) 

2,1 
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Consequently there exists a nonnegative summable sequence {cj{t) = Cj{u, v, i))jgz whose 
sum is 1 such that for all j G 

l|A,-(^;.Vn)||^i^. < 2-^-(^-i) [\,{t) l^uo\\u\\j^s+r+C'^\\vf,, \\u\\^s-A dr. 

Jo y ' ^2.1 ' / 

(3.35) 

Similarly we obtain thanks to Proposition 13: 



\v.Vq\\^s-i <c\\v\\ d \\q\\ss^^ ^'^ll'^ll.^ + + 



and there exists a nonnegative summable sequence whose sum is 1, once again denoted 
by {cj(t) = Cj{q,v,t))j^^, such that for all j € Z 



+ C^^me.x{l,^)\\v\\\ (IklUf-i+^lkb^J )^^^- (3-36) 

^ ^21 J 

Writing the Bony decomposition, 

\\v.Vq\\ss^^ < \\f,Vq\\^s^^ + \\fygV\\^s^^^ + \\Riv,Vq)\\j;,s^^ 

Using (4.125), we get that: 



\\Tvqv\\Bs^ + \\Rivyq)\\Bs^ < c\\v\\ ,^m\bs^ < cibll (IklU-i 

^2,1 ' ^2,1 ^ ' 

which implies that there exists a nonnegative summable sequence whose sum is 1, once 
again denoted by (cj(t))jgz such that: 

u2^\\A,{fy,v + R{v,VqmLiL2 < C2-^'(^-i) f Cj{T)\\v\\ . (\\q\\^si+u\\q\\^ )dT. 

JQ ^2 1 ^ ' 

(3.37) 

Unfortunately, we are not able to estimate this way T^Vq: there are too many derivatives 
involved in high frequency for the density fluctuation. But if we restrict to the low 
frequencies when j < we can deal with this term: as the frequencies of Sk-iV-A^Vq 
are localized in an annulus of size 2'^, there exists an integer Ni (only depending on the 
parameters chosen in the Littlewood-Paley theory) such that: 



and for all j < 0, 



Ajit^Vq) = ^j{Sk-ivAkVq) 

\k-j\<Ni 



2^'\\Aj{f„Vq)\\L2 < 2^' Yl lbllL-2'=||Afcg||i2 

\k-j\<Ni 
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There exists a nonnegative summable sequence {di{t) = di{q,t))i^i- , with < 1 

and a constant C = C(iVi) such that for ah I < Ni, 

where the low frequency summation is well-defined because: 



^ /<Ari Y l<Ni Y l<Ni 

As we consider j < 0, we have j — Ni < k < j + Ni < Ni < — log2 e if e is small enough, 
so that we can write 



and then 

where d'j{t) = ^\k-j\<Ni dk'2!^^~^'^^ is in as a convolution of summable sequences 

(and its norm is bounded and only depends on s). Finally, we obtain that for all j < 0, 

< 2-^-(-i)^*4(r) (^J^u^q\\^s^.,s + ^ll^lli^^llgll^.J dr (3.38) 
For all J G Z and t £ I, let us introduce: 

Uj{t) = \\A,u\\L^L^+uo2'^\\A,uhiL^ + {l+u2^) (^\\Ajq\\L^L^ + umm{^,2''^)\\A,q\\LiL?j 

(3.39) 

and 

u{t) = hllz-B-i+lkllzrBri'^'''''^''^rBia+''°''''''^t^S'+''''^ 

(3.40) 

d 

Using estimates (3.35), (3.36), (3.37), (3.38) in (3.33) (and the fact that ^ L°°, see 
appendix), we obtain that there exists a nonnegative summable sequence whose sum is 
1, denoted by {cj{t))j^z such that for all j < 0, 

U,{t) <C^ C/,(0) + (1 + u2^)\\AjF\\l^l, + ||A,G||^ii2 

' 

+ C^:|2-^-(-i) rc^.(^)[(n,ax(l,i3) + -)IIKT)||% +C\\v{T)\\,d^Au{T)dT . 

(3.41) 
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When J > we are not able to estimate the problematic term Ty'Vq with this method. 
The key idea is to get rid of every advection term using a Lagrangian change of variable 
as in [24] to [27], [15] and [6]. This is the main difficulty of the article and the object of 
the following section. 

3.2 Lagrangian change of coordinates 

As stated before, the aim of this part is to get rid of the advection terms involved in 
system (LR^). Let us first consider the localized equations (as usual we set fj = Ajf...) 
written in the following way: 

dtQj + Sj^iv.Vqj + div tij = fj, 

dtUj + Sj-iv.Vuj - Auj + pVqj - ^V((/>e * qj - qj) = gj, 

where the external force terms are defined by: 

fj = Fj + (^Sj-iv.Vqj — Aj^v.Vq)^ and gj = Gj + (^Sj-iv.Vuj — Aj^v.Vu)^ . 

These two terms can be estimated thanks to the following commutator estimate from 
[15] (we refer to lemma B.l from appendix B): 

Lemma 1 ([15]) There exists a sequence {cj)j^i G /^(Z) such that ||c||;i(2) = 1 and a 
constant C = C{d, a) such that for all j G Z, 

II J i J n niL - 1 II "s^^nL- ^2.1 



Remark 18 Let us emphasize that we considered, up to adding external force terms 
that we can control (fj and gj), the advection by the low frequencies of v, that is exactly 
the quantities that we were not able to deal with. 

Let us set ipj^t as the flow associated to Sj-iv: 

(dt^Jj^tix) = Sj-iv{t,'ipj^t{x)) 

we can also write: 

ft 



V'j,i(x)=X+ / Sj^lv{T,'4jj^r{x))dT. 







Thanks to propositions 14 and 15 from the appendix (we refer to [15] or [6]), there exists 
a constant C such that: 



llff o V'i.i IIlp ^ e II^IIlp for all function g in , 



l-D'^Vtilli- < C2('=-i)j(e^^ - 1) for A; > 2, 
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def 



where 

V{t)"= I \\Vv{T)\\LoodT. (3.44) 
Jo 

Let us also emphasize that in the present case, the jacobian determinant of the change 
of variable will play a crucial role in the obtention of uniform estimates with respect to 
£ (contrary to the case of lemma 2.6 from [3] where it produces a term that we are not 
able to sum): 

fdei(L'Vi,t(a;)) = ei'o('^'"^^-i^)("'^^-.-(-))'^", 

[detiDll^-lix)) = e-/o(div5,-it;)(r,X,(r,M))dr ^ g- /o*(div 5,_i«){r,^,>o,^-l(x))dr ^ (3-^^) 

where Xj (r, t, x)) denotes the two parameter flow associated to Sj-iv (we refer to (4.126) 
in the appendix). 

We will now perform the announced change of variable, for a function h, let us define 
h = ho ipj^t = h(t,7pj^t)- Then we have dtqj{t,x) = {dtQj + Sj-iv.S/qj){t,7pj^t{x)), which 
provides the following system: 

dtQj + divuj = fj + R], 

^ K ^ ^ ^ n o (3.46) 

dtUj - Auj + pVqj - -^V{(j)e * Qj - Qj) = Qj + Rj + Rj + kRj, 

where most of the remainder terms Rg, R^ and Rg are exactly the same as in [6] (with 
the same convention: if / : M*^ — t- M™" is a differentiable function then Df denotes the 
Jacobian matrix of /, and V/ is the transposed matrix of Df.): 

R]{t,x) := Tr{Vuj{t,x) • (/^ - V^-,^(^j- t(x)))) , 
R]it,x) := Vqj{t,x) ■ {VijjJ{^j,t{x)) - Id) 

and i?| := fiR'j + (A + ^f)i?| with 

Rf\t,x) := Tr((V^-,nV'j,t(x)) - h) ■ VDu%x) ■ D^l^-j{i;j,t{x)) 

+VDu]{t, x) ■ {Dijjjiijj^tix)) - Id)) + Vuiit, x) ■ AV^-/(^,- t(x)) 
Rf\t,x) := Tr(^Duj{t,x) ■ OkDi^jJU^j^tix))) 

+ Yl ^bcu]{i, x) ■ dkij'^'^ipj^tix)) ■ daip~l'\'ipj,t{x)) 

a,b,c,b^a,cj^k 

d 

+ Ydlfi]{t,x) ■ (^dkipjl'''{t,'il;j^tix)) - Id) ■ di^jj'\t,'iljj^t{x)) + {di'4)jy{ipj^t{x)) - h)). 

i=l 

There is only one additionnal remainder term compared to [6]: 

R, = ° - (3.47) 

Thanks to the definition of (hr and 6 we obtain that: 



*/-/ _! 

£2 £2 



{z)ifix-£z)- fix))dz (3.48) 
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and the difficulty comes from the fact that we need to obtain bounds in LjL^ that are 
uniform with respect to e, that go to zero when t is small (estimated by e'-'^^*) — 1), 
and that do not involve too many derivatives (that is 2'^-^) for the density fluctuation 
in high frequencies. Remember that formally "^^i^zJ! _i. go we formally have 

— )• Aqj o ipj^f — A(qj o V'j^t)"- this term has no reason to be small in e and the best 
we can hope is to get uniform bounds with respect to e. 

Dealing with Rj is the object of the rest of this section. 
3.3 Precisions on the capillary term 

Before dealing with Rj let us briefly go back in this section to the convolution term 
written in (3.48): for a function /, 

Proposition 7 For any suitable function f and any s G M, the following two norms are 
equivalent: 

^s+2,s ~ II 11^,^^ . (3.49) 




Proof: When we consider a frequency localization of /: if j G 



iO = -2 A,/(0, 



£2 £2 

and thanks to the Plancherel formula, 

Thanks again to the fact that on M+, g : x ^ l — e~^ is increasing and h : x ^ {l — e~^)/x 
is decreasing, if cq2^ < < Cq2^ , we can write: 

< < 

and 



£2 £2 



Cl £2 - £2 - ^2 £2 

therefore 

Remark 19 Obviously we can replace cq by any positive constant. 

On the other hand, it will be useful to compare ^— ^ — 5 and min(^, 2^^). 
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If 2^ > \ then > Cq and as function g is increasing and bounded from above 

by 1: 



1 - e-^o 1 - e-^o^'^^^ 1 . 1 



22j 



< 5 < ^ = min( — , 22J). (3.50) 



If 2-' < ^ then 2-^e^Co < Cq and as function h is decreasing as for all x > 0, 
< h[x) < 1 we have: 



cl - e^cp^^ - 



and 



1 —c'^e 2''^ 1 

(1 - e-^o)22i < < 012"^ = d min(^, 22J). (3.51) 



From (3.50) and (3.51), we deduce that for all j G Z: 

(1 - e-"o)min(^,22j) < „ < max(l, eg) min(^, 2^^), 

and 



min(^,22^-)||Az/||^2. 



Multiplying by 2-'^ and summing over j £ Z, we obtain that 

,,4*/-/,, 



II ^2 llB|_i ll-/ llsr^'"- " 

Remark 20 This new formulation seems more natural than (1.1) as instead of a Exed 
frequency threshold there is a continuous transition zone between the parabolically reg- 
ularized frequencies and the damped frequencies. 

Remark 21 Notice that the norm obtained by replacing <f>{^) = e"'^'^ by e""'^'^ is 
equivalent (the multiplicative constants depending on a). 



Remark 22 Our space is a particular case of the general hybrid Besov spaces introduced 
by R. Danchin (see [11, 3]): we refer to the appendix for the fact that ijl"*"^'* ~ Be 



2 

e 



Remark 23 In the U-setting, we can prove that for all j G Z: 

„ 0,*A,/-A,/ ^^^^ ^ Cmin(l,22^-)||A,/||,. 
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3.4 Estimates on the capillary term 

In this section we wish to focus on the capillary term: 



= ^ ° V-i,* ^2 

Thanks to 

Vqj = V{qj o ipj^t) = yqj o ^j,t X Dipj^t = ^qj o ^j,t x {Dipj^t - h) + Vgj o V'j,*, 
we obtain the following decomposition: Rj = Ij + Ilj with 

Ij = * ^2 — ~ where gj = Vqj o ^j^^ x (/^ - -DVj,*), 

^ (3.52) 
_ * Vgj - Vqj _ (pe * (Vgj o ^j,t) - Vgj ° i^j,t 

^~ £2 °W ^2 

The main difficulty of this paper is then clearly established and consists in estimating 
(locally in frequency) the commutator between the Lagrangian change of variable and 
the non-local operator: 

LeU) = ^1^4—^ = 1 / <l,{z) (fix - ez) - fix)) dz. 

For a function /, and for all j G Z we set fj = Ajf and: 

If, = Il',{f) = 'tllJyll o _ <t>e*{f,oi.,^-f,oi.^, _ ^^^^^^ 

Theorem 5 Let cr G M. There exists a constant C = Ca^d such that for all f € Be^'^''^ , 
there exists a summable positive sequence {cj{f))j^z whose sum is 1 such that for all t 
so small that 

e2^^-l<^. (3.54) 

and for all j £ TL, 

||//;.(/)k. < Ce^^(y + e^^^ - l)c,(/)2-^-|| '^- * {" , 

£" 2,1 

where Vit) = \\\/v{T)\\L°°dT. 

Remark 24 As a by-product we obtain that under the previous assumptions, if t is 
small enough, 



-||(</>,*A,/)oV^,,,-0,.(A,/oV>,^i)||^2 <Ce^^(y + e^^^-l)c,(/)2-^- f^ [ ^ ||^.^ 



even if neither of the left-hand side terms are spectrally localized. 
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Proof: Let us first rewrite the non-local operator: 



Le{f){x) = '^'*{ h x) = 4 / 4>e{x - y) ifiy) - fix)) dy. (3.55) 

In the works of T. Hmidi, S. Keraani, H. Abidi, M. Zerguine (we refer to [24], [25], [26] 
and [27]) the key idea is to express the difference II j as an integral formulation and 
try to retrieve the desired Besov norm thanks to an equivalent expression of this norm 
involving finite differences of / of order 1 (that is expressions of the type T-yf — f where 
T-yf{x) = f{x + y)) or order 2. 

Remark 25 Due to the equivalent forms of our hybrid norm, it would be useless here 
to simplify the expression of Hjif) and write IIj = (((/>£ * fj) o ^^^^ — (/j o i^j,t))- 

We refer to the following results of [3] for these useful alternative characterizations of 
Besov norms (we also refer to proposition 1.37 with a simpler proof in the specific case 
of Sobolev spaces): 

Theorem 6 ([3], 2.36j Let s g]0, 1[ and p,r £ [l,oo]. There exists a constant C such 
that for any u £ C'f^, 

n-^w II ^ II W'^-y^ - ^Wlp II ^^11 II 



'Ivl" 



and in the case where s = 1, 



Theorem 7 ([3], 2.37) Let p,r £ [l,oo]. There exists a constant C such that for any 



u£C'^, 



r'-ill II ^ II W'^-y'^ + '^y^ - '^Alp h n 

II^IIbi,, - II [^] llL'-(Rd;y^) - ^ll^llsi,; 



Remark 26 We emphasize that in the second case, when the regularity index s is an 
integer, we have to use finite differences of order 2 instead of order 1. In the present 
article, due to the capillary term we will have to completely rewrite these results. 

For example, in [25], the authors need to estimate in L^ the commutator |Z)|°(/ o tp) — 
{\D\°' f)oip in terms of the 5^ ^^-norm of / (with a s]0, 1[), where ij) is the flow associated 
to a divergence- free vectorfield. To do so they write both terms in the commutator in 
a unified shape and manage to get the result thanks to a non-local expression of the 
fractionnal derivative: 

Jm.<>- \x — y\ ^ 

and thanks to an equivalent expression of the Besov norm, using finite differences of 
order 1. Let us emphasize that in this case, the vectorfield is divergence- free therefore 
the Jacobian determinant of ip is 1, which simplifies the estimates. 

In our case, we have to estimate J/j which is a commutator between the non-local 
operator and the Lagrangian change of variable. Our operator is still non-local but the 
vectorfield is not incompressible anymore, and more important, we will have to construct 
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an analogous of the finite difference expression of our hybrid Besov norm. As the Jacobian 
determinant of ipj^t is not constant we cannot use the following estimate, used in [24], 
[25], which is due to Vishik (see [37] in the divergence- free case): 

\\Al{Ajao^p)\\LP < C2^l-''-^l || VV''^^"^-'-') • \\Aja\\LP (3.56) 

In the non measure-preserving case we can recall the following result: 

Lemma 2 ([3] lemma 2.6) There exists a constant C > such that for all global diffeo- 
morphism of M'^ i/j, for all p > 1, for all tempered distribution a, and for all j,l G Z, we 
have 

\\Al{AJao^p)\\LP < C2-J'|| JV^-^IUco ||A,-a||Lp(||VJV'-'||L- || JV'IIl- + 2'||V^||loo). 



Using this result in our non measure-preserving case involves an additionnal term (due 
to the Jacobian determinant) that we are not able to deal with (indeed, the first term 
in our case is not summable in /). On the other hand, we cannot hope to use a better 
estimate than the following one, that is small when t is small but consumes too many 
derivatives: 

\\fMA-))-fj(.-)\\L^<cv{t)2^\\f,\\L. 

To bypass this difhculty we will simply write directly the expression in the Lagrangian 
variable and precisely trace the Jacobian determinant. Let us go back to our estimate: 

^^■(^) = if / Mi^jA^) - y) ifjiy) - fM^A^))) dy 

- [ - y) UA^jAy)) - fA^jA^))) dy) ■ (3.57) 

Let us now study II'Ai^jJ{x)) instead of //j(x): 

ii'Ai^-Jix)) = i f / - y) if Ay) - f'ji^)) dy 

U^jlix) - y) ifA^jAy)) - fA^)) dy) ■ (3.58) 
Performing in the second integral the change of variable y = 'iIjJI{z), we obtain: 

/ Mi^-Ji^) - y) ifA^iAy)) - fA^)) dy 

= [ Ui^j}{x)-i^j}{z)){fAz)-fA^))\detDijj}{z)\dz (3.59) 

Thanks to (3.45) we have (recall that Xj is the two-parameter flow associated to Sj-iv: 
Xj{T,t,z) = ipj^roip-liz)): 

t 

(div Sj^iv){T, Xj{T, t, z))dT 



det{D^j}{z)) = e Jo 
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so that 



<Pe{x - y) {fj{y) - fj{x)) . 



(^eii>-l{x)-i^^l{y)) - (divS',„i?;)(r,Xj(r,t,y))dr 
1 — — e -^0 

i'eix - y) 



dy). (3.60) 



Recall that for all x G W^, <Pe{x) = j2^^& ^ and perform the change of variable 
z = 2—^, we obtain: 



1 



\y\ 



e 4 (/.(^_ey)_/.(:c)). 



2 / |V>-i(x)-V^i(x-ey)p~ 



1 - e 



e .^0 



(div S'j_iu)(T, Xj(r, t, X — ey))dT 



dy 



(3.61) 



As we want to estimate the norm of this quantity, that is a Besov norm with integer 
regularity index s = 0, the finite difference of order 1 will not be sufficient for our need, 
and we will have to introduce finite differences of order 2. Indeed, using the present 
quantity would only involve a term in e\y\ and when estimating in low frequencies, there 
would be either a multiplicative coefficient 1/e or an additionnal derivative term (that 
would prevent any convergence when —j is large). To be able to do a correct estimate 
we need at least e^|yp. 

To do this we simply write IIj = + and perform the change of variable 

z = —y in the second integral. If we set: 



A : 

B 
C 
D 



\ylf_ \i^-l{x)-^J^{x-Ey)f 
4 I e^lyp 



(div5j_iw)(r, Xj(T, t,x — ey))di 



(3.62) 



4 \ e^lyp 



(div Sj^iv){T, Xj{T, t,x + ey))dT, 
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then 



+ / e 4 (/^.(x + ey) - fj{x)) [1 - e^e^] 



Illj + IVj. (3.63) 



where 



(3.64) 

Remark 27 Even if there is only a finite difference of degree 1 in the second term, we 
wiU he able to get correct estimates thanks to the coefRcient e^e^ — e'-'e^ . 

We will make an extensive use of the following elementary consequence of the mean- value 
theorem: 

Lemma 3 For any x,y gM., — e^| < \x — y\e^^^^^'y\ 
Let us begin with IIIj: 

Proposition 8 Under the previous assumptions, there exist a positive constant C = C^^^d 
and a nonnegative sequence [cj = Cj{f))j£z whose summation is 1, such that if t is so 
small that e^*^^(*) — 1 < |, we have: 

WIIIjllL^ <C{V + e -l)e 2 ■> Cj\\ 



To prove this result we will successively prove the following lemmas: 

Lemma 4 There exists a constant C such that for all j € Z, /, and all t is so small that 

2CV(t) _ l< 1 

xl/ e--^\\fj{.-ey) + f,i. + ey)-2f,{.)\\L2dy, (3.65) 



where V is defined in (3.44). 
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Lemma 5 For all a £ M there exists a constant C^^d such that for any f G -Bf i' t-^sre 
exists a nonnegative summable sequence {cj{f))j£z with ||cj(/)||;i(g) = 1, such that 

4 / e--^\\f,{.-ey) + f,{. + ey)-2f,{.)\\L2dy 

< C^4'^-'''c,{f)\ &*i~^ \\B^, ■ (3-66) 



Proof of lemma 4: we only need to estimate with lemma 3 the following decomposition 
of the coefficient appearing in Illj : 

l — e e =1 — e +e (1 — e j. 

The first estimate is straightforward: as \B\ < C Jq HVS'j-if (r)||Loo(fr < CV{t), 

|1 - e^l < \B\e\^\ < Cy(t)e^^(*). (3.67) 

Before estimating the second term, as in the work of T. Hmidi and S. Keraani ([25]) we 
need to give a precise bound for the variation ratio of the flow: 

Lemma 6 Under the same assumptions, for all x, y G M'^ with y ^ 0, we have: 



< e2^^(*) - 1. 



Proof: For ah x / y G R'^, thanks to (3.43), 

\x — y\ ^' 

Similarly, for all x 7^ y G M*^, 

|x-y| 

and applying it at il)~l{x) and '4}jl{y), 

|^-i(x)-V'-i(y)| - 
so that we get that for all x and y 7^ 0, 



e2|y|2 



in other words, 



1 < _,-2Cl/W < 1 ( \i^ji{x)-i^j}{x-ey)\'' 
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that is exactly the announced result. ■ 

Coming back to lemma 4 and thanks to lemma 3, we immediately get that: 



Gathering (3.67) and (3.68), we get 



/ IM^('^,2C\/(t) _ i \ I |2 



so that, taking the L^-norm in the expression of II Ij from (3.64), we obtain: 

Wllljhl < ^,^\, 2 / e 4 ||/^.(._ey) + /^.(. + £y)_2/,(.)||^.||l-e^e^||i^dy, 

(^C7y(t)e"^ + M!e-¥e^(^^""'*^-^)(e2^^W - 1)) 

X ||/,(. - ey) + /i(. + ey) - 2f,{.)\\L2dy. (3.69) 



1 



Assume that t is so small that e^*^^^*^ — 1 < ^5 then (using again that for all x > 0, 

xe~^ < |e~^/^) 

— e — i-e~^'' * < liLLg — 8-<-e-T^, 
4 4 e 

which ends the proof of lemma 4. ■ 

Remark 28 Similarly, writing L^{fj) = ^{L^{fj) + L^{fj)) and performing the change 
of variable z = —y in the second term lead to: 

\y\^ 

^^(/i) = \d^2 / ^ ^ - ^y) + hi^ + ey) - 2/i(a;)) c^y, 

taking the L^-norm we obtain that: 

_\yl 

\ &*%~^' \\L-<:rj:^ [ e A \\f^{. - ey) + f,{. + ey) - 2fj{.)\ydy, 
and as Tai^jf) = Aj^Taf) (which is localized in frequency), we immediately get that: 

_\yl 

\ &\{~^ \\b^,<^ I e lQ\\f{.-ey) + f{. + ey)-2f{.)\\ dy. 
On the other hand, lemma 5 gives the reverse estimate so that: 

_\yl 

||g||^.+2,. ~1 / e 16\\f{.-ey) + f{. + ey)-2f{.)\\^.dy. 
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Proof of lemma 5: let us adapt the proof of [3] (th 2.36 and 2.37). As fj = Ajf, we 
can write that for all x and y, 



fjix - ey) + f,{x + ey) - 2fj{x) = [r.^yA^f + T,yA,f - 2Aj/j (x) 

and as in the proof of theorem 2.37 from [3], thanks twice to the mean value theorem, 
we get that: 

\\T^ey^,f + T,yA,f-2Ajf\\L2<C2^^e^\y\^ J2 W^ffh^, 

li'-il<i 

and thanks to the definition of the hybrid Besov norm || ^^-p— ^ ||^<t = || / 11^^+2, <t, there 
exists a nonnegative summable sequence (whose summation is 1) {cj)j£z such that: 

||r_.A-/ + -.A-/-2A,/||,.<C2^V|,|^ E W^l^n 'l^^''^.-.- 

\j'-j\<i ^'P' > 

Moreover 

1 1 max(e2,2-2i') 



min(^, 22i') min(^, 22i) max(e2, 2-2i) ' 
and thanks to the following result, 

Lemma 7 ([11], Proposition 5.3j Let a > 0, a, 6 G M. Then we have 

max(a,2-") ^ fl if a>b, 
max(a, 2"'') " |2^-" if a<b. 

as \j' — j\ < 1, we obtain that: 



mm(^,22J') - mm(4,,22J)' 

so there exists a nonnegative summable sequence of summation 1, that we will also denote 
by {cj)j^z, such that: 

||r_,,A,/ + T,yA,f - 2A,f\\L2 < C„2^^e^\y\^2-^''cj max{e^ 2-^n\ & * \ 
On the other hand, we can also roughly write: 



iT^ey^jf + TeyAjf - 2A,-/||i2 < 3||A,/||i2 < C2-^- cj max{e\2-^^ )\\^^^^^^\\ . 



"2,1 

Combining the last two estimates, if we define j^y £ ^ the unique integer satisfying 
1 < ^\y\2^^^ < 2 (that is e|y|2-'=« ~ 1) , we get the following statement: 

IKsy^jf + T,yAjf - 2Ajf\\L2 < C^2~i''cj max(e2, 2-2i)|| '^^ ^^'^ bf,, 

f22..2|,|. if,|,|2.<l (,<,,^) 
ll ife|2/|2J>l {3>3ey). 
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We can now prove lemma 5: let us follow the proof from [3] and begin by splitting the 
integral: 

\ I e--^\\f,i.-ey) + f,{. + sy)-2f,{.)\\L2dy 

e \ Je\y\2J<l Je\y\2i>l J 

X C,2-J% max(e2,2-2J)||^^^^^^||^^. (3.71) 



We have: 



1 



/ \y\ e 16 dy = C / r^e isr'' dr < Cd / re 32dr = Cd(l-e 32!^) 
J|y|<4T "'0 Jo 



<C,min(l,-2^). (3.72) 



Similarly, 

/■ |j^|2 /■OO 2 /"OO 2 X 

/ e'^dy = C I e'Tar'^-^dr < Ca re'^dr = Cde~^^^^ . (3.73) 
Plugging this into (3.71), we obtain that: 



\ I e--4\\fji.-ey) + f,{. + ey)-2f,i.)\\L2dy 



< a2'^-c, max(l, -i^)|| '^^ * ^ ||^. (2^%^ min(l, -l^) + e-^ 

< C^,d2-^"c,-|| '^"*{~^ ||^. (i?i + iZa) , (3.74) 

£: 2,1 

where 

r 1 1 

= 2^%^ min(l, ^2^) max(l, ^2^) = 1, 

_ 1 1 _ 1 _ 1 

i?2 = max(e sil^W,^— iSI^W) < max(l, Ce e^J^W) < C, 

which ends the proof of lemma 5. ■ 

Lemmas 4 and 5 immediately imply proposition 8, and we can now turn to IV f 

Proposition 9 Under the same assumptions, there exist a positive constant C = Co-,^ 
and a nonnegative summable sequence {cj = Cj{f))j^i whose sum is 1, such that if t is 
so smaU that e^'^^^*) — 1 < 5, we have; 

\\IVj\\L^ <C^^d[V +e -l)e 2 J CjW . 



To prove this result we will prove the following lemmas: 
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Lemma 8 There exists a constant C such that for all j G /, and all t so small that 
2CV{t) _ i< 1 

mm{l,e2^\y\)e-^-^\\f,{. + ey)- f,{.)\\L2dy, (3.75) 

where V is defined in (3.44). 

Lemma 9 For all a £ M there exists a constant Ca^d such that for any f € B21, there 
exists a nonnegative summable sequence (cj(/))jgz with ||cj(/)||;i(2) = 1, such that 



4 / mm{l,e2^\y\)e-^\\fj{. + ey) - f,{.)\\L2dy 



<a,,2-^%(/)||^i^^^||. (3.76) 



Proof of lemma 8: as announced here we can only rely on a finite difference of order 1, 
and we need to carefully estimate the following coefficient to obtain summable sequences 
in low frequencies: 

e^e^ - e^e^ = e^(e^ - e^) + e^(e^ - e^). 

Let us begin with the second term: thanks to lemma 3 (we refer to (3.62) for the expres- 
sions of the various terms involved), 

- e^l < e--(l^l'l^l)|i) - B\ < (I^M^I)(|Z)| + \B\) 

<2e^y [ \\div Sj^i{T)\\L^dT <2Ve^y. (3.77) 
Jo 

and then, 

\yL(f>2cv[t) _ 

e^\e^-e^\<2e 4 'Ve^^. (3.78) 

This estimate is useful only for high frequencies: indeed for the low frequency regime, 
as in the proof of proposition 8, after integration with respect to y, the result won't be 
summable when j goes to —00. In order to get a suitable estimate for low frequencies, 
let us use once again the elementary mean value theorem: 



J ^div Sj-iv{T, Xj{T, t,x + ey)) — div Sj-iv{T, Xj(T, t,x — ey))^ dr 

^^cv{t) J 2e\y\\\V (divSj-iv{T,XjiT,t,.))^\\L^dT. (3.79) 
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Obviously: 

V (div Sj-iv{T, Xj{T,t, z))j = Vdiv Sj^iv{T,Xj{T,t,z)) x DXj{T,t,z)), 
so that, thanks to the Bernstein lemma (we refer to the appendix), (3.43) and (3.45): 

||v(div5,_i^;(T,X,(r,t,.))) < \\Vdw Sj^Mt, Xj{r,t, .))\\L^\\DX,{r,t, 

< ||Vdiv5j_i?;(r)||Looe^^W < 2^\VviT)\\L^e^^^'\ (3.80) 

and ^ 

|e^ - e^l < e'^^^2e2^y\ [ \\Vv{T)\\L^dT < e^^^ 2e2^ \y\V (t) . (3.81) 

Jo 

Combining (3.78) and (3.81) we end up with the better estimate: 



\yl 

e^|e^ -e^l < Ce 4 e^^V{t)mm{l,e2^\y\), (3.82) 



which now can be used for all frequency. The same has to be done for the other term: 
rough estimates first imply that 

\y^( 2CV{t) _ -^^| |2 

|e^-e^|<e--^(l^l'l^l)2max(|^|,|C|)<2e 4 \yL^e^cV(t) _ (3_83) 



And more precisely (needed later for the low frequencies convergence), we write that: 

\e^-e^\ <e 4 ^ (3.84) 

with 

\yl ( \i^-h^)-i^J,h^ + ^yT _ \i^jjix)-i^lhx-ey)\' \ 
4 1^ e2|y|2 ,2\y\2 j 

= 4^2 (iV',"'!^) - ^Jji^ + - IV'-'(^) - i^-Jix - ey)\^) . (3.85) 

Thanks to the identity |ap — |6p = (a + b\a — b), where (.|.) is the usual scalar product 
in M'^, we get: 

1^1 " 4^2 {^^7}^^^ ~ ^7,t^^ + ^2/) - i'iiix - ey) ipjjix - ey) - ipjjix + ey)^ 
- 4^2 1 2V',"' (x) - i^jj {x + ey) - ^jj {x - ey)\ x \ iPjj (x - ey) - ^PjJ (x + ey)\. (3.86) 

Thanks to the mean value theorem (used twice for the first factor and once for the 
second), we can write: 

\F\ < ^(,|y|)2||Z)2^-i||^^ X ie\y\)\\Dij-lU^ < \e\yf\\D%^X-mj,l\\L-, 
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and using the estimates for the flow (see (3.43)), we obtain that 

|F| < Ce\yf2^e^^{e^^ - 1), 

and ^ 

I^L Cp2cy(t) _ IN 

|eA _ <ce ^e|y|32^e^^(e^^ - 1). (3.87) 

From (3.83) and (3.87) we deduce that: 

\y^( 2cv(t) _ 

eS|e^-e^|<Ce^^W(e^^W-l)min(l,e2^|y|)|ype 4 . (3.88) 

Finally, thanks to (3.88) and (3.82), 

lle^e^-e^e^llijo <Ce^^W(y(t) + e^^W-l)min(l,e2^|y|)(f + |2/|2)e 4 >. 

(3.89) 

Taking the L^-norm in (3.64), and thanks to the previous estimate, we have: 
^/ e 4 min(l,e2^|2/|)(l + |y|2)e 4 ^ ^ ||/^.(. + ey) _ 

(3.90) 



and as previously, if t is small enough so that 6^*^^^*) — 1 < ^, 

x^f + 8 min(l,e2^|y|)||/,(. + 6y)-/,(.)||i2dy (3.91) 

We conclude thanks to the fact that for all x > 0, xe~^ < -e~2. ■ 

— ' — e 

The proof of lemma 9 will follow the lines of the proof of theorem 2.36 from [3]: 



fj{x + ey) - fj{x) = (^r,yA,f - A,/j (x) 

and thanks to the mean value theorem, 

\\T,yA,f-Ajf\\L.<C2^e\y\ J] ||A,,/||^2, 

li'-il<i 
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and as in the proof of lemma 5, using the definition of the hybrid Besov norm, there 
exists a nonnegative summable sequence (whose summation is 1) once more denoted by 
(cj)jgz (depending on / and t) such that: 

\\TsyAjf-A,f\\L2 < C72^e|y|2-^''^c,-max(e2,2-2^')||^^^^^||^.^. 
On the other hand, a rough estimate provides: 

\\TeyA,f - A,f\\L2 < 2\\A,f\\L2 < C2-^-c,-max(£2,2-2j')||^i^4— ^IIb? ■ 
so that we end up with: 



\TsyA,f - Ajf\\L2 < a2-^'%-max(e2,2-2^')||^^^^-^||. 



\2^yl its|,|2^<l (;<;.„) 
\l ifE|9|2J>l {j >],,}. 



Returning to the integral, we have: 
1 



£2 



min(l, 82^ \y\)e-^ ||/,(. + ey) - fji-^L^dy 

<max(l,^)a,,2'^-c,(/)|| '^^ * ^ \\^.^ J^^mm{l,e^2'^\y\^)e-'4dy. (3.93) 
Then the same computation as in the proof of lemma 5 leads to: 
1 / mm{l,e2^\y\)e~^\\f,{. + ey) - fj{.)\\L2dy 

<C^,,2-^'^c,{f)\ & *i~^ \\B?,- (3-94) 

This concludes the proof of lemma 9 and together with lemma 8 it implies proposition 9. 
Finally, propositions 8 and 9, with the first estimate from (3.43) end the proof of theorem 
5. ■ 

3.5 End of the proof 

Once obtained the desired estimates of the additionnal remainder term introduced by 
the effect of the Lagrangian change of coordinates on the non-local capillary term, we 
can go back to the use of the linear estimates from proposition 2 on system (3.46): for 
all / G Z, as qj{0) = ^^(0, o(.)) = qj{0), 

W'^iUjUl^l^ + 1^0 2^' 1 1 A; 11^1^2 + (1 + 1^2^) ^11 A/gj- 11^00^2 + i/min(^, 2^')||Ajgj||iii2^ 
(1 + u2^)\\Aiq,{0)h2 + ||A,n,(0)||i2 + (1 + u2^)\\Aif, + AiR]\\^l, 



<C ^2 



+ WAiQj + AiR) + AiR) + KAiRjWi^iL^ 



. (3.95) 
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Let us recall that thanks to (3.41), all we need is to estimate the high frequencies, that 
is {qj,Uj) for j > 0. For this, as in [15] and [6], let us define some Nq G Z (that will be 
fixed later), and write: 

\ l>j-No 

Moreover, as in [15] and [6], we can use the version of (3.56) given by lemma 2 in the 
non measure-preserving case (we refer for example to [3] lemma 2.6, or [15], lemma A.l) 
to obtain: 

WSj^NoQjWL^ < Ce^^ (e^^ - 1 + 2^^0e^^) \\qjh2, 

so that: 

lkillL2 <Ce^^ [(e^^-l + 2-^0e^^)||g,-|U2+ ^ \\A1qjh2] , (3.96) 
V l>j-No J 

and for the velocity, 

\\uj\\L2 < Ce^^ I (e^^ - 1 + 2"^«e^^)||7Xj||i2 + W^i'^ih^ ] ■ (3-97) 

\ l>j-No ) 

Multiplying (3.96) by 2-' and using that in the summation, I > j — Nq, we obtain: 

2^lk,||L2<Ce^M(e^^-l + 2-^o,cv^)2.||,^.||^,+ 2'2^o||A,^,||,. ) . 

\ l>j-No J 

Then, going back to Uj (we refer to (3.39) for the definition), we can write that for all 
J >0, 



Uj{t) < Ce 



cv 



[e^v _ 1 _^ 2-^°e^^) U,{t) + max(l, 22^«) ^ 



l>j-No 



KiUj\\L^L2 + VQ2^^\\l\iUj\\^L2 + (l + i/2')(^||Azgj||i^^i2 + i/min(^,22')||Ajgj||ii^2)^ 

(3. 

where we used once again lemma 7 to prove that for all I > j — Nq: 

min(l,22^) = min(l,2^') "^^"|f r,7, < 2^^o min(l, 2^^ (3.99) 
e'^ max(e"', 2~"'-?) 
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Thanks to (3.95) we obtain that for all j > 0, 



Uj{t) < Ce 



cv 



1 + Ujit) + max(l, 2^^'')C ^ x 



((1 + ^2')\\Aiq,{0)h2 + \\Aiu,iO)\\L2 + (1 + (||A;/,-||^i^. + || A^i?] 



l>j^No 



L2J 



, (3.100) 



Except Rj, all of the remainder terms are the same as those from [6]. As we will sum 
for I > j — Nq, we need to provide estimates involving some positive power of 2^^K To 
do this we use the well-known method of Vishik (see [37] ) which takes advantage of the 
Bernstein lemma: the idea is to derivate the function in order to obtain a positive power 
of 2^~^ which is summable over I > j — Nq. Let us detail this for example on fj: thanks 
to (3.43), 

\\Ajj\\L2 < C2-^\\AiVfj\\L2 < C2-'||V/jO V,-,j||L2||Z?V'i,i||L» 



and thanks to lemma 1 (for o" = s — 1), we can write that there exists a nonnegative 
sequence {cj{t))j^i G such that ||c||;i(2) = 1 and a constant still denoted by C = 

C{d, s) such that for all j € Z, 



\Aif,\\L2 < C2^'-'e^^ + c,(t)2-^'(^-^)||V 



V\\ d 



■n-J -"2,1 

Similarly, we obtain that (using lemma 1 for o" = s in the second case): 



(3.101) 



\Aigj\\L2 < C2^-'e^^ WGjWl^ + Cj{t)2-^('-^^\\Vv\\ a \\u\\^s-i 



2'\\Aifj\\L2 < C\\AiVfj\\L2 < C2^-'e^^ 2^\\Fj\\l2 + Cj{t)2-^^'-^^\\Vv\\ , \ 

\ ^2,1 



With the same argument we obtain (we refer to [6], section 3.2 for details): 



IBS 



(3.102) 



\\AiR]\\l2 < C2^~\e^'^ - l)e^^2muj\\L2 < C2J-'(e^^ - l)e^^22i||n^.||^2 
21\\AiR]\\l2 < C\\AiVR]\\l2 < CT~\e^^ - l)e^^2''^\u,\\L2 
\\AiR]\\l2 < C2^--He^^ - l)e^^2^||(?,||i2 
\\AiR';\\l2 < C{\\ + Ail + ^i)2^-\e^^ - l)e^^22i||n,||i2 



(3.103) 



Remark 29 JVote that the last estimate from the first line is valid only when j > 0. Let 
us recall that we already treated the low frequencies thanks to the linear estimates (we 
refer to (3.41) ). 



43 



All that remains is to estimate kRj: as for the other terms, a direct estimate will not 
be sufficient as we eventually have to sum over I > j — Nq. We then need to use the 
Bernstein lemma in order to be able to factor by a positive power of 2^~K Let us recall 
that Rj = Ij + IljiVq) with (we refer to (3.52)): 

= — — |i — ^ where gj = Vqj o ^j^t x {h - Dipj^t), 

Jj'Jf) = 't>^*fj-fj o J. 0e*(/j°V'j,t)-/joV'j,t 

Let us recall that section 3.4 is devoted to the following estimate (see theorem 5 for 
details): there exists a nonnegative summable sequence Cj{t) (whose sum is 1) such that 
for all j G Z (recall that V{t) = \\Vv{T)\\L-odT), 

mmiL^ < CMe^^(y + 6^^^ - l)c,{t)2"n &*{-^ \\^^ , 
Let us begin with Ij: a rough estimate gives: 

ll^/^jllL2 = II ^lll,2 < Cmm(-2,2 )||5illL2 

< C min( 1 ) 1 1 V(7j o V^,. f 1 1 ^2 1 1 (/^ - D V',. < ) 1 1 . (3 . 104) 

Thanks again to lemma 7, 

min(l,2^-) ^ min(1.2^^)""°";;'^-;;' < min(l. 2'^^) ( \„ , " " ^ ' £ J'- 

V2' ' ^£2' Vax(e2,2-2') - ^2' ^|22('-J) ifj</, 

so that (thanks to (3.43)) as max(l, 22('-j)) < 2^^02^(1-3) [Nq will be taken positive), 
||A,/,||^. < Ce^^ie^^ - l)22^022('-^-)min(l,22^-)||Vg,||^., 

which implies that in our case, in order to use Vishik's trick we need to derivate three 
times A//j because we have to absorb 22('~-') and obtain a positive power of 2-'"' (which 
is summable over I > j — Nq): 

We then compute the third derivatives of gj = Vqj o ipj^f x (Id — Dipj^t) and obtain that: 

l|V^5,llL2<Ce^^^(e^^-l)2=^^||V,,||^., 
which allows to conclude that 

\\AiI,\\l2 < C7e^^^(e^^ - l)2^'-'22^° min(l, 2^^)11 Vg,- 11^2. (3.105) 
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For the second term (//j(V(7)), we have to derivate one time and decompose Vl/j into 
three parts: 

V//j = o iPj^t ■ D^j^t — 

= Aj + Bj +Cj, (3.106) 

where 

= ^ ° ^j,t ■ {Di)j,t - Id), 

Bj = Il'^iV^q), 

^ _ (t>e * (v^gj o i'j,t ■ Wj^t - Id)) - v^gj o i^i,t ■ {D'^j,t - Id) 

and thanks to the Bernstein lemma: 
\\AiIIj\\L2 < C2-'||A;V//,-||i2 < C2-' (\\AiA,\\l2 + \\AiB,\\l2 + \\AiC,\\l2^ . (3.107) 

Thanks to theorem 5 with a = s — 2 and for t small enough (e^*^^ — 1 < ^, see the proof 
of theorem 5) there exists a nonnegative summable sequence (cj(r))jgz whose sum is 1 
such that 

HA R II ^ IIR II <r (M ^J'CV (yj2„ ^\o-j(s-2) \\ 't>£ * " V^g n 

\\A.iBj\\j^2 < \\Bj\\]^2 < Ls,de [V+e - i)Cj[\/ q,T)2 -^"^ '\\ ^ lls^-^i 

< C,,,e^^(y + e2^^ - l)c,(v2g,r)2^-2-^-(^-i)||^i^^|^||^.-i. (3.108) 
A direct estimate gives: 

||A;^j||i2 < \\Aj\\L2 < ll^^^^^-^l ° V'i,t||L2 • WD^pj^t - Id\\L-° 

< Ce^^ie^^ - 1)2^ II '^^ * ^% ~ ^'^^ \\l2, (3.109) 

and finally, for the same reason as in Ij, rough estimates of AiCj will provide min(p-, 2^') < 
22{'-i) min(^, 2^-') and we once again have to derivate Cj: this time we will derivate only 
two more times because we already have 2~' in factor: 

2-^||AiCj||i2 < C2-^^\\AiV^Cj\\L2 < Ce^^^ie^^ - 1)22^«2J-' min(-^, 22j)||Vgj||i2. 

Plugging this together with (3.109) and (3.108) into (3.107), and using (3.105) allows to 
get the estimate on Rj: 

WAiRjW < Ce^^(V + e2^^ - l)2^--'2^^°c,(t)2'^-(-'i)|| * " ||^.-x (3.110) 

S 2,1 
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Plugging (3.101), (3.102), (3.103) and (3.110) into (3.100), then summing over I > j - Nq 
implies that for all j > and for t small enough (e^*^^ — 1 < ^) we have: 



Uj{t) < Ce 



2CV 



^^cv _ 1 + 2-A.o) U,{t) + 25^oC ^ (1 + ^2^')|k,-(0)||^. + ||^,(0)|U. 



+(lW)||F,-||^i^. + ||G,-||^i^. + (e^^-l) ((l + |A + H+^ + z^)2'-'h,||ii^. + ||Vc7,-|Lii2) 



+ / 2-^-(-i)c,(r)||Vz;(r)|L. (IkWb.-i + Hk(r)||^. ^ + ||n(T)||^.-i)dr 



Bo 



2,1 

+ 4(y + e2^^-l) f c,{r)2-^^^-')M-llllzll,. 

'0 



. (3.111) 



As j > 0, we have 1 < min(i, 2^^) and 



(1 + |A + /.I + + ^)22^||n,|Li^. + llVg.ll^.^. < ( ^ + \^ + + ^ + J_ | u,it), 



so that (remember that t is so small that e*^^ — 1 < A' 



C/,(t) < -C 



(e^^ - 1) f 1 + 2^^oc7 ( ^ + 1^ + /-'+^ + - + J,)) + 2-^o) jj,[t)+ 
25^oC ^ [ (1 + u2^)\\qmh2 + ||u,(0)|U. + (1 + ^2^')||F,||,i^. + ||G,|| m^. 



+ / 2-^-(-i)c,(r)||V^;(r)|| (|k(r)||^. i + ^||g(r)||^ + ||n(r)|| i)dr 

JO -02,1 



+ 4(y + e^^^-l) / c,(r)2-^-(-i)..2|l^l^^i^ll^ 

JO 



I/^ ./n £ 2,1 



(3.112) 



Let us fix A'^o > large enough, and take t small enough so that: 

■-C-2-^«<i, 
4 - 8 

^C(e^^ - 1) f 1 + 2^^^C .2 ( ^ + " + ^1+^ + " + ^)]<\ 
Then we obtain that for all j > 0, 



(3.113) 



u^{t) < 3C25^°c ^ (1 + u2^)\\q,mL^ + iwrnh^ 



+ (1 + u2^)\\Fj\\l, , + \\Gj\\l1l2+ / 2'^^^-'hj{T)\\Vv{T)\\ , U{T)dT 

Jo 



2CV 



1) / c,(r)2 
^0 



-i(^-i)j.2||^ 



* Vq — Vg 



IB; 



.-idr]. (3.114) 
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Now, if t is so small that: 



3C2^^'>C ^{V + e^^^ -1) <l- — , (3.115) 



then for all j > 0, 



2^^'-^^Uj{t) < 3C2^^°C ^ ( 2J'(^"i)c/j(0) + (l + i/2^')2^'(^^i)||Fj||iii2+2J'(^-i)||Gj||iii2 

+ fc,iT)\\VviT)\\.. [/(r)(ir ) + /■^^.(^)ii^llZpZ^II d^. (3.116) 
JO -82,1 / ^ JO £ ^'^ 

Going back to the low frequencies, if we take K = (2C ^2 in (3.41), then for all j < 



2^-(-i)C/,-(t) <C^ 2^-(-i)C/,(0) + (1 + z.2^-)2^-(-i)||A,F||^i^2 + 2^-(^-i) || A,G||^j^2 
+ 2-^-(-i) rc,(T) ((max(l,i3) + l)||^;(r)f + ||^(r)|| ] [/(r)dr 

+ i2--''(^-i) / Cj(r) (zyo||n||^s+i dr. (3.117) 

2 Jo V 2,1 e £ / 

Summing over j € Z gives that 



+ ( ^ + " + ^'+^ + % max(l,J3)) /V'(r)f/(r)^ (3.118) 
1^0 '^•^ Jo / 



where 



V{t) r ||VT;(r)||i^dr < W{t) A||VKt)|| , + HrW , )rfr. (3.119) 



|2 ^ 

^2^1 ^2^1 



and thanks to the Gronwall lemma, we obtain that for t small enough (satisfying condi- 
tions (3.54), (3.113) and (3.115)), 



U{t) < 2C ^ (c/(0) + ||F|| .15=- + HI^IIlJbi , + \\G\\ 



2C „,(l±lA±i4±ii±i: + „,ax(l.l))H'(t) ^ ^ 
X e ^0 (3.120) 

Using conditions (3.54), (3.113) and (3.115) gives the result for small times t (it does 
not depend on the initial data, only on V{t) and such that (3.113) is satisfied), then we 
globalize the result (for example as in [26], end of section 4) for all time t by subdividing 

[0, t] into intervals [Tj+i, Tj] where we have e^^ ■^'^^ V{T)dT _^ < i as well as the analogous 
of (3.113) (we also refer to [3] chapter 10 for a connexity argument). ■ 
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3.6 Extension of the results 



In this short section, we mention the following extension of the linear estimates from 
proposition 2 in Besov spaces constructed on U spaces with r 7^ 2. This high-frequency 
result can be obtained from the Fourier expression of the solution of the linear system as 
done in [6] (see section 2.1). 

Proposition 10 Under the same assumptions as in proposition 6 (adapted to r ^ 2), 
there exists a constant such that for aU j > jg, 



Using the Lagrangian change of variable, remark 23 and the methods from [6], we 
obtain the analogous in the L''-setting of theorem 1.5 for high frequencies, which is similar 
to the key proposition 6 from [6] that leads to the main result of [6]. Let us recall that 
in this case, the regularity index of the Besov spaces is ^ — 1 which is negative when 
d < r, that allows initial data with large modulus provided that they have fast enough 
oscillations. 

4 Appendix 

The first part is devoted to a quick presentation of the Littlewood-Paley theory and 
specific properties for hybrid Besov norms used in this paper. The second section to 
general considerations on flows. 

4.1 Besov spaces 

4.1.1 Littlewood-Paley theory 

As usual, the Fourier transform of u with respect to the space variable will be denoted by 
J-{u) or u. In this section we will briefly state (as in [6]) classical definitions and properties 
concerning the homogeneous dyadic decomposition with respect to the Fourier variable. 
We will recall some classical results and we refer to [3] (Chapter 2) for proofs (and more 
general properties). 

To build the Littlewood-Paley decomposition, we need to fix a smooth radial function 
X supported in (for example) the ball i?(0, |), equal to 1 in a neighborhood of i?(0, |) and 
such that r 1— )• x('^-er) is nonincreasing over ]R_|_. So that if we define ip{S,) = x(?/2) — x(0' 
then (f is compactly supported in the annulus € IR'^jCq = f < |C| ^ Co = f } and we 




(3.121) 



have that 




(4.122) 
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Then we can define the dyadic blocks {Ai)ii=z by A; := ip{2 ^D) (that is A^n = (p[2 ^Cju{^)) 
so that, formally, we have 

u = Y^Kiu (4.123) 
I 

As (4.122) is satisfied for ^ 7^ 0, the previous formal equality holds true for tempered 
distributions modulo polynomials. A way to avoid working modulo polynomials is to 
consider the set 5^ of tempered distributions u such that 

lim llS^nllioo = 0, 
>— 00 

where Si stands for the low frequency cut-off defined by Si := x(2~'D). If n G 5^, (4.123) 



is true and we can write that Siu = AgU. We can now define the homogeneous Besov 

k<l-l 

spaces used in this article: 

Definition 2 For s € M and 1 < p, r < 00, we set 

\u\\^s '■= ( 2'"^^|| A^nll^p^ if r < CO and := sup 2'''||A/ii||LP. 

We then define the space B^.^ as the subset of distributions u & S'^ such that \\u\\^s 
finite. 



IS 



Once more, we refer to [3] (chapter 2) for properties of the inhomogeneous and homoge- 
neous Besov spaces. Among these properties, let us mention: 

• for any p G [1, 00] we have the following chain of continuous embeddings: 



p,oo ' 



d 



• if p < 00 then B^^ is an algebra continuously embedded in the set of continuous 
functions decaying to at infinity, in particular we make in this paper an extensive 

d 

use of the injection B21 L°°; 

• for any smooth homogeneous of degree m function F on M'^\{0} the operator F{D) 
maps Bp J. in Bp~/^. This implies that the gradient operator maps 5^ ,, in Bp~r^. 

We refer to [3] (lemma 2.1) for the following result describing how derivatives act on 
spectrally localized functions: 

Lemma 10 (Bernstein lemma) Let < r < i2. A constant C exists so that, for any 
nonnegative integer k, any couple {p,q) in [l,oo]^ witli q > p > 1 and any function u 
of L^, we have for aiJ A > 0, 

Suppu c B{0,XR) =^ \\D''u\\l', < C*-'+^A''+^^^"i^||u||Lp; 
SuppSc{eGM^/rA< 1^1 < i?A} ^ C^'^-^A'^llnllLP < \\D^u\\lp < C'^^ X'' \\u\\ lp . 
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This implies the fohowing embedding result: 

Proposition 11 For all s £M, 1 < pi < P2 ^ oo and 1 < ri < r2 < oo, the space Bp^ j.^ 

.s-d(-L_J_) 

is continuously embedded in the space Bp2^r2^^ ''^ 



In this paper, we mainly work with functions or distributions depending on both the 
time variable t and the space variable x. We denote by C{I;X) the set of continuous 
functions on / with values in X. For p S [1, oo], the notation L^{I; X) stands for the set 
of measurable functions on / with values in X such that 1 1— )• belongs to LP{I). 

In the case where / = [0, T], the space L^([0, T];X) (resp. C{[0,T]; X)) will also be 
denoted by L^X (resp. CtX). Finally, if / = we alternately use the notation L^X. 

The Littlewood-Paley decomposition enables us to work with spectrally localized 
(hence smooth) functions rather than with rough objects. We naturally obtain bounds for 
each dyadic block in spaces of type Lj^L^. Going from those type of bounds to estimates 
in Lf^Bp j. requires to perform a summation in ^""(Z). When doing so however, we do 
not bound the Lf^Bp ^. norm for the time integration has been performed before the 
summation. This leads to the following notation (after J.-Y. Chemin and N. Lerner in 
[9]): 

Definition 3 For T > 0, s G M and 1 < r, p < oo, we set 

One can then define the space L^^Bp ^. as the set of tempered distributions u over (0, T) x 
R*^ such that lim„^_oo S„u = in LP([0,T]; L°°(R'^)) and \\u\\rp < oo. The letter T 

is omitted for functions defined over M"''. The spaces Lf^Bp^. may be compared with the 
spaces Lj,Bp^^ through the Minkowski inequality: we have 

ll'ullrpRs <||'u||rPRs \i r>p and ||n||rpr}s >\\u\\jpp,s if r < p. 

All the properties of continuity for the product and composition which are true in Besov 
spaces remain true in the above spaces. The time exponent just behaves according to 
Holder's inequality. 

Let us now recall a few nonlinear estimates in Besov spaces. Formally, any product 
of two distributions u and v may be decomposed into 

uv = TuV + T^u + i?(n, v), where (4.124) 

TuV := ""^^ Si-iuAiv, TyU := ""^^ Si^ivAiu and ii(ti,f):=^^ AiuAi'V. 
I I I |/'-«|<i 

The above operator T is called "paraproduct" whereas R is called "remainder". The 
decomposition (4.124) has been introduced by J.-M. Bony in [5]. 

In this article we will frequently use the following estimates (we refer to [3] section 2.6, 
[11], [18] for general statements, more properties of continuity for the paraproduct and 
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remainder operators, sometimes adapted to Lf^B^ j. spaces): under the same assumptions 
there exists a constant C > such that: 

II^Ii'^IIrs < C||u|ktx) IIuIUs < C||n|| d IbllRs , (4.125) 

-°2,1 ■°2,l r"? -°2,1 

^2,1 

\\tuv\\^s+t < C||u||^^ ^ ^ < C7||n|| , ^ < 0), 

^2,1 

\\R{u,v)\\^s,+s2 <C\\u\\^n J\v\\^s2^ <C||n|| d||i;||^^2^ (si + S2>0), 

-°2,1 

\\R{u,v)\\ _d<C\\R{ 

'^T^)\\^^l+'>2 < C'll^'ll^^i,^ ||^'||^^2^ (si + S2 > 0). 

4.1.2 Complements for hybrid Besov spaces 

As explained, in the compressible Navier-Stokes system, the density fluctuation has two 
distinct behaviours in some low and high frequencies, separated by a frequency threshold. 
This leads to the definition of the hybrid Besov spaces. Let us begin with the spaces that 
are introduced by R. Danchin in [11] or [3] (we will use these spaces only in the appendix 
to prove estimates with the Hybrid norms introduced in (1.1)): 

Definition 4 For q > 0, r € [0, oo] and s € M we denote 

\\u\\gs,r =^ ^ 2'* max(a, 2"')^"^ ||Ain||i2 
i&z 

For example with r G {l,oo}: 

\\u\\§s,oo = ^ 2'('*~^)||A,n||i2 + ^ a2'''||A,ii||i2, and 

Klog2{^) '>log2(^) 

II^IIb-^= E 2'(^"''^||A;n||i2+ Yl ^2'^l|A/u|U2, 

^<log2{^) i>log2(^) 

Remark 30 As stated in [3] we have the equivalence 

— (\\u\\ns~l +a\\u\\r)a 1 < ||u||gs,oo < ||n||r>s-l +a||u||r>s . 

We refer to (1.1) and Proposition 7 for the precise expression of the hybrid norm used in 
the present article. Let us just mention that this particuliar hybrid norm is accurate for 
our problem, but it is also related to the hybrid norms introduced by R. Danchin: 

||.||^.+2,. = ||.||^,,2 

Let us now state the following result (Proposition 5 from [7]) 
Proposition 12 Let s G M, a > 0. For all q G Ba°^ n Ba^ , we have 

\\q\\%s ^ Ikll 5"-°° Ikll R^'i 
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Remark 31 For all q G B^J^ n ^| ^ = Bl'°° we have 

\\q\\§s.oo < + Ikllsi 

and when e > is small enough, for all q € Be~^^'^, we have 

so we can use the hybrid norms introduced in (1.1) and we will in fact use the following 
results: 

Proposition 13 Let s € M. There exists a constant C > such that for all < e < 1, 

and aU q G B^^ n S| ^ n Bl+^^'~^ n 51+^'", we have 

\\q\\ls^ < c(lkb-i + \\q\\bsJ{\\q\\bi+','-^ + 



4.2 Estimates for the flow of a smooth vector-field 

In this section, we recall classical estimates for the flow of a smooth vector-field with 
bounded spatial derivatives. We refer to [15] or [6] for more details. We also refer to [24] 
for the incompressible Navier-Stokes case. 

Proposition 14 Let v be a smooth globally Lipschitz time dependent vector-field. Let 
W{t) := Jo \\Vv{t')\\L^ dt'. Let ipt satisfy 

iptix) = x+ [ v{t',^t'{x))dt'. 
Jo 

Then for all t G M, the flow ipt is a smooth diffeomorphism over and one has ift>0, 

||I?Vt^lL=o < e^W f\\DMt')\\L^e'^(''^dt', 
Jo 

WD'^t'h- < e^(*) 1^ \\DMt')\\L^e''^^''Ut' + S^e^W \\DMt')\\L-e'^'-''^dt'y . 

As in [6] we also introduce the two-parameter flow (t, t' ,x) i— t- X{t, t' ,x) which is (uniquely) 
defined by 

X{t,t',x) = x+ [ v{t",X{t",t',x)) dt". (4.126) 
Jt' 

Uniqueness for Ordinary Differential Equations entails that 

X{t, t",X{t",t', x)) = X{t, t',x). 
Hence Vt = -^^(^,0, •) and ip^^ = X{0,t,-)- 
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Proposition 15 Under the previous notations, the jacohian determinant of X satisfies: 

det{DX{t,t',x)) = e/t'(div^')(r,X(r/,x))dr^ ^^-^27) 

and 

det{D'4>t{x)) = e/o(div«){r,i/.,(x))dr^ 

Proof: differentiating (4.126) with respect to x, one gets by virtue of tlie chain rule, 

DX{t,t',x) = Id+ [ Dv{T,X{T,t',x)) ■ DX{T,t',x)dT. (4.128) 
Jt' 

This immediately implies that: 

dt{DX){t,t',x) = Dv{t,X{t,t',x)) ■ DX{t,t',x), 

and 

dtdet{DX{t,t',x)) = tr {Dv{t,X{t,t' ,x))) ■ det{DX{t,t' , x)), 
so that we obtain the result. ■. 

The authors wish to thank Raphael Danchin, Taoufik Hmidi, Miguel Rodrigues and 
the anonymous referees for useful remarks and suggestions. 
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